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Existence of the thermodynamic limit and asymptotic behavior of 
some irreversible quantum dynamical systems 

Abstract 

This dissertation discusses the properties of two open quantum systems with a 
general class of irreversible quantum dynamics. First we study Lieb- Robinson bounds 
in a quantum lattice systems. This bound gives an estimate for the speed of growth 
of the support of an evolved local observable up to an exponentially small error. In 
a second model we study the properties of a leaking cavity pumped by a random 
atomic beam. 

We begin by describing quantum systems on an infinite lattice with associated 
finite or infinite dimensional Hilbert space. The generator of the dynamics of this sys- 
tem is of the Lindblad-Kossakowski type and consists of two parts: the Hamiltonian 
interactions and the dissipative terms. We allow both of them to be time-dependent. 
This generator satisfies some suitable decay condition in space. We show that the dy- 
namics with a such generator on a finite system is a well-defined quantum dynamics 
in a sense of a norm-continuous cocycle of unit preserving completely positive maps. 

Lieb- Robinson bounds for irreversible dynamics were first considered in [10] in the 
classical context and in |20] for a class of quantum lattice systems with finite-range 
interactions. We extend those results by proving a Lieb-Robinson bound for lattice 
models with a more general class of quantum dynamics. 

Then we use Lieb-Robinson bounds for a finite lattice systems to prove the exis- 
tence of the thermodynamic limit of the dynamics. We show that in a strong limit 
there exits a strongly continuous cocycle of unit preserving completely positive maps. 

V 



Which means that the dynamics exists in an infinite system, where Lieb-Robinson 
bounds also holds. 

In the second part of the dissertation we consider a system that consists of a beam 
of two-level atoms that pass one by one through the microwave cavity. The atoms 
are randomly excited and there is exactly one atom present in the cavity at any given 
moment. We consider both the ideal and leaky cavity and study the time asymptotic 
behavior of the state of the cavity. 

We show that the number of photons increases indefinitely in the case of the 
ideal cavity. In case of leaky cavity, we prove that the mean photon number in the 
cavity stabilizes in time. The limiting state of the cavity in this case exists and 
it is independent of the initial state. We calculate the characteristic functional of 
this non-quasi-free non-equilibrium state on the Weyl algebra. We also calculate the 
energy fiux in both the ideal and open cavity and the entropy production for the 
ideal cavity. 
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Chapter 1 
Introduction 

1.1 Introduction 

The study of open quantum systems is significant to various areas of application 
of quantum theory such as quantum optics, quantum information theory, atomic 
physics and condensed-matter physics. In this dissertation we assume an open system 
dynamics of a quantum system. First we study Lieb-Robinson bounds for a lattice 
system with irreversible time-dependent dynamics. Then we study the properties of 
a state of the microwave cavity pumped by a random atomic beam. The dynamics 
of the system is also time-dependent and irreversible. 

In the past years Lieb-Robinson bounds have been shown to be a powerful tool to 
turn the locality properties of physical systems into useful mathematical estimates. 
Lieb-Robinson bounds provide an estimate for the speed of propagation of signals in 
a spatially extended system. 

Lieb-Robinson bounds were so far considered on a lattice system (for example 
Briefly, if A and B are two observables supported in disjoint regions X and F, 
respectively, (see section 12.2.21 for the definition of the support) on a lattice system 
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and if Tt denotes the unitary time evolution of the system, a Lieb-Robinson bound is 
an estimate of the form 



\\[Tt{A),B]\\ < C7e-^('^(^'^)-"l*l\ (1.1.1) 

where C, fi and v are positive constants and d{X, Y) denotes the distance between X 
and Y. Here v is called Lieb-Robinson velocity. See Chapter 12.2.71 for more details. 

From this inequality we see that for the small times |t| < d{X, Y)/v the norm on 
the right-hand side is exponentially small. The reason we consider the commutator 
on the left-hand side of the Lieb-Robinson bounds is the following. The commutator 
between the observables A and B is zero if their supports are disjoint. The converse 
is also true: if the observable A is such that its commutator with any observable B 
supported outside some set X is zero, then A has a support inside the set X. This 
statement is also approximately true in the sense: suppose that there exist some e > 
such that II [A, i?] II < e||i?|| for any observable B that is supported outside a set X 
and some observable A. Then there exists an observable A^"^^ with the support inside 
the set X that approximate the observable A: \\A — A^'^'^W < e. 

So Lieb-Robinson bounds say that the time evolution of the observable A with the 
support in a set X is mainly supported in the 5-neighborhood of X, where 6 > v\t\ 
with V is the Lieb-Robinson velocity. 

In quantum information if Alice has access to the observable A with support X 
and sends a signal to Bob who has access to the observable B with support Y then 
the strong enough signal that Bob can detect propagates with the Lieb-Robinson 
velocity v. 

In 1972 Lieb and Robinson first showed the existence of the light-cone such that 
the amount of information signaled beyond it decays exponentially [12]. This result. 
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known as Lieb- Robinson bounds, has been improved in [9], [10], [H] -[18] , for exam- 
ple, making the bound on the speed independed of the dimension of the spin spaces 
and extending the class of Hamiltonians describing the system. In these works the 
Lieb-Robinson bounds was considered for the unitary evolution, when the dynamics 
is described by the Heisenberg equation. 

We consider a more general situation, when the dynamics is described by a Marko- 
vian dynamical semigroup (see Chapter 12.2. 4l for more details) with a time dependent 
generators. Lieb-Robinson bounds for irreversible dynamics were, to our knowledge, 
first considered in [10] in the classical context and in [20] for a class of quantum 
lattice systems with finite-range interactions. We extended these results by proving 
the Lieb-Robinson bound for a long-range exponentially decaying time dependent 
interactions. 

One of the main applications of Lieb-Robinson bounds is the existence of the 
dynamics in the thermodynamic limit. The existence of the thermodynamics limit 
is important as a fundamental property of any model meant to describe properties 
of bulk matter. In particular, such properties should be essentially independent of 
the size of the system which, of course, in any experimental setup will be finite. Our 
results are described in Section [3l 

In the second part of the thesis we use the results about the irreversible dynamics 
to study the properties of a leaking cavity pumped by a random atomic beam. 

The micromaser system that we consider consists of a beam of two-level atoms 
that pass one by one through the microwave cavity. The atoms are excited with some 
probability p. The cavity is modeled by a single photon mode. The beam is tuned 
in such a way that there is exactly one atom in the cavity at any given time. While 
in the cavity the corresponding single atom is able to interact with a single mode of 
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the cavity field. 

The interaction between the atom and the cavity in our model leaves the state 
of the atom invariant which may be interpreted as the limit of heavy atoms and soft 
photons. Hence, we focus our study on the time asymptotic behavior of the cavity 
state. 

We do this both in the case of a perfect (non-leaky) cavity as well as in the case 
of a leaking cavity modeled by adding a Lindblad-Kossakowski dissipative term to 
the generator of the Hamiltonian dynamics (see section 12.2.41 for details). We show 
that in the case of the perfect cavity, the expected number of photons in the cavity 
increases indefinitely in time if and only if < p < 1. In other words, fluctuations of 
the atom state are necessary for the pumping to occur. We derive an expression for 
the mean photon number for any time. 

When the leakage is taken into the account, the mean- value of photons number 
stabilizes in time. The limiting state of the cavity exists and it is independent of 
the initial state. We study this non-quasi-free non-equilibrium limiting state by 
computing its characteristic function (see section [2.1.3l for the definition of quasi- free 
and non-equilibrium steady state). We also calculate the energy flux in both the ideal 
and open cavity and the entropy production for the ideal cavity. 
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1.2.1 Lieb- Robinson Bounds and the Existence of the Ther- 
modynamic Limit for the Class of Irreversible Quan- 
tum Dynamics 

In Chapter [3] we prove Lieb- Robinson bounds for a class of the irreversible quantum 
dynamics and show the existence of the thermodynamic limit of the dynamics. 

We consider a general situation, when the dynamics is described by a Markovian 
dynamical semigroup. The generator of the dynamics consists of both Hamiltonian 
and dissipative terms that may depend on time. In the finite volume A the generator 
L is of the following form 



for all local observables A. Here for Z C A, Z) is the local time dependent inter- 
action describing Hamiltonian terms and Z) are local time dependent bounded 
operators. We put an exponential decay condition in space on the The detailed 

set up is discussed in section ISTTl 

Fix T > and, for all observables A with finite support in A, let v4(t),t G [0,T] 
be a solution of the initial value problem 




(1.2.1) 



ZcA 



^z{t){A) 




d 

di 



A{t) = Umit), AiO)=A. 



(1.2.2) 
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For < s < t < T, define the family of maps {7f^s}o<s<t C B{Aa,Aa) by ■y^s{A) = 
A(t), where A{t) is the unique solution of fll.2.2p for t G [s,T] with initial condition 
A{s) = A. In section 13.21 we show that the dynamics 7^ generated by (11.2. ip exists 
and forms a norm-continuous cocycle of unit preserving completely positive maps. 
In section 13.31 we prove the Lieb- Robinson bound of the form 

\\^iHB)\\ < C(/C, B)e-M'^{^'^)-Ii|)^ (1.2.3) 

where /C is a completely bounded linear operator on the space of observables with 
support X that vanishes on 11 (see section 13.11 for the definition of the completely 
bounded map), B is an observable with support Y, d{X,Y) is the distance between 
the supports X and Y, C{IC,B) is constant that depends on /C and B and v is 
a Lieb-Robinson velocity. It is important that the operator /C can be taken in a 
form }C{B) = [A, B], which makes the original Lieb-Robinson bound for reversible 
dynamics (ll.l.ip a particular case of the Lieb-Robinson bound in the general form 

In section 13.41 we look at the dynamics 7^ when the set A becomes infinitely 
large. We prove that the thermodynamic limit of the dynamics exists in the sense 
of strongly continuous one-parameter cocylce of completely positive unit preserving 
maps. 

1.2.2 Non-equilibrium state of a leaking photon cavity pumped 
by a random atomic beam 

The system we consider consists of the one-mode microwave cavity and a beam of 
randomly excited atoms. The cavity is described by quantum harmonic oscillator 
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with the Hamiltonian 

He = eb*b, 

where e > and b*, b are boson creation and annihilation operators (see section [2.1.11 
for more details on these operators). 

The beam of atoms is described as a quantum spin chain 

Ha = Y1 EVn, 

n>l 

where E > and for any n > 1, the atomic operator rjn '■= {cr^ + ll)/2, where is 
the third Pauli matrix. 

Atoms are randomly excited with some probability p. They enter the cavity 
successively and there is only one atom present in the cavity at any given time. 
Atoms in the excited state interact with the cavity and produce a quantum field. 
The interaction between the n— th atom and the cavity is the following 

K{t) = X[{n-l)r,nr]{t){X{b* + 6) ® Vn), 

here Xi{^) is a characteristic function of a set /, which is here to make sure that 
the n—th atom interacts with the cavity only when it is present there, which happen 
during the time interval [(n — l)r, nr]. The detailed set up of the system is given in 
the section urn 

At first we assume that the cavity insulates an atom-photon system from de- 
cohering interactions with its environment. So the Hamiltonian for the system is the 
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sum of the Hamiltonian of the cavity and atoms and the interaction between them 



H{t) = eb*b ® 11 + ^ 11 ® + ^ 

X[(n—l)T,nT] 

n>l n>l 

In section 1121 we show that in this case the expected number of photons in the cavity 
oscillates around a linearly increasing in time function and only a beam of randomly 
exited atoms (0 < p < 1) is able to produce a pumping of the cavity by photons. 

In a more general situation we should not assume that the cavity is perfectly 
insulated. Photons may either escape the cavity or be absorbed into the cavity walls 
at some constant non-zero rate a > 0. 

In this case the evolution of the system is described by a Markovian dynamics 
with the following generator 

mips) = -t[H{t),ps]+cxb{ps)b* - ^{b*b,ps}, 

where ps is a state of the system and cr > 0. Here for any observables A and B, 
{A,B} := AB + BA denotes the anticommutator. 

In section l4l3] we show that the expected number of photons in the cavity stabilizes 
in time, which means that the energy leaking out of the cavity equals to the energy 
coming into the cavity. 

In section 14.41 we compute the characteristic function of the limiting state of the 
cavity on the Weyl algebra (see Section [2. 1.2p . The limiting state is an infinite convex 
combination of quasi-free states, and, in general, is not quasi-free itself (see Section 
12.1.31 for the definition of a quasi-free state). The limiting state of the system is a 
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non-equilibrium state in a sense that it is not in a form 
for any /3 G C. 

Section 14.51 is devoted to the calculation of the energy flow for the perfect cavity 
(14.5. ip and the leaking cavity fl4.5.2p . The entropy production for the ideal cavity is 
presented in Section fl4.5.3p . 
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Chapter 2 
Preliminaries 

2.1 Quantum systems 
2.1.1 Fock space 

The quantum-mechanical states of a particle form a complex Hilbert space "H, for 
example it could be L'^{R'^). The Hilbert space of n identical but distinguisable 
particles is the tensor product 

7/" = (g)7{. 

i=l 

If the number of particles is not fixed, the states are described by vectors in the Fock 
space J^{T-L) which is given by 

^(H) = (g)7/^ 

n>0 

where "H^ = C. Therefore a vector ijj G J-'i/H) is a sequence of vectors tp = }n>0) 
where ijj^"'' = fi®...®fn^ G "H for 1 < < n. In quantum physics identical particles 
are indistinguishable which is reflected by the symmetry under the interchange of the 
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particle coordinates. 

The first case is when the components ip^"'^ of each ip are symmetric under the 
interchange of coordinates. Such particles are called bosons, they form the Bose-Fock 
space J-'^{'H). The second case is when the components ip^'^'^ of each ip are anti- 
symmetric under interchange of each pair of coordinates. Such particles are called 
fermions and they form the Fermi-Fock space J-'_('H). 

Bose- and Fermi-Fock spaces are defined with the help of the operators P+ and 
P_ which are defined on J^{'H) as follows, for all /i, .../„ G "H 

P+(/i ® ... ® /„) = 5^ A, ® ... ® 

TT 

where is the sign on the permutation n. The sum is taken over all permutations 
TT : {1,2,.. .,n) — (vTi, 7r2, ...7r„). These operator are extended by linearity to two 
densely defined operators with ||P±|| = 1. Then extending P± by continuity one gets 
two bounded operators of norm one. The Bose- and Fermi-Fock spaces are given by 

A number operator is defined on J^{T-C) by 

= {nV;(")}„>o 

with the domain D{N) = {ip = {^(")}„>o : En>o^^ll^^"^ll^ < 

The particle creation and annihilation operators on J^{'H) are defined as follows. 



2.1. Quantum systems 

For each f e H, a{f)tlj^^'^ = and a*(/)^/;(°) = / and 



12 



aAf)ifl ® - ® fn) = V^(/, /l)/2 ® ... ® fn, 

a>(/)(/i ® ... ® /„) = + 1/ ® /i ® ... ® /n, 

for any fkE'H,k&Z. Extending by linearity one gets two densely defined operators. 
If ^/'("^ G we get 

||a^(/)ril < V^ll/lllini, l|a>(/)|| < v/;^||/||||V^(")||. 

Therefore ajr(f) and a*jr{f) have well-defined extensions to the domain of D{N^^'^) of 
Ari/2. 

From the definition one can see that a* operator is indeed the adjoint of a. For 

all 0, ^ G D{N^/'^) and f en we have that 

(a>(/)0,V^) = (0,a^(/)7A). 

The creation and annihilation operators on the Fock spaces J^±{'H) are defined by 

a±U) = P±a{f)P± = aif)P±, 
alif) = P±a*if)P± = P±a*if). 

In future when it is clear what Fock space we are working in we may drop the 
indices ±. In Section H] we will use one particle boson creation and annihilation 
operators b*, b that act on the Hilbert space "H. 

The important properties of creation and annihilation operators are canonical 
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commutation relations (CCR) for boson creation and annihilation operators and 
canonical anti-commutation relations (CAR) for fermion creation and annihilation 
operators. 

Boson operators satisfy the following CCR relations 

[aM),^+i9)] = 0=[a*{f),a*{g)], (2.1.1) 
[a+(/),a*(^)] = (/,<?) 11. 

Fermion operators satisfy the following CAR relations 

{a4f),a4g)} = = {a*_(/), a^^?)} (2.1.2) 
{a4f),a*_{g)} = {f,g)l. 

2.1.2 Weyl operators 

In this section we will consider the Boson- Fock space J-'+('H). So we will drop the 
plus index on the creation, annihilation operators. 

It is convenient to introduce the family of operators {b{f); f G H} as follows 

b{f) = a{f) + a*{f), 

where a*, a are boson creation and annihilation operators. Each b{f) is a self-adjoint 
linear operator on the Fock space. The canonical commutation relations (12.1.11) in 
terms of the operators b now take form 



[b{f),b{g)] = 2zlm{f,g). 
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W{f) = exp{zb{f)}. 

The operators b{f) as well as the creation and annihilation operators are un- 
bounded, while the Weyl operator is unitary. This is one of the reasons it is easier to 
work with Weyl operators than with creation and annihilation operators. 

The Baker- Campbell-Hausdorff formula shows that for two operators X and Y 



if [X, [X, Y]] = = [Y, [X, Y]]. Using this formula one can get the Weyl form of the 
canonical commutation relations 

W{f)Wig) = e-5i"^(^'5)vr(/ + g)= e-'^'^^^^'^^Wig)W{f). 

In Section H] the creation and annihilation operators are b* and b that act on the 
Hilbert space T-L. The Weyl operator is then defined as follows: for every a E C, 



W{a) = exp{ab-ab*}. (2.1.3) 
From Baker-Cambell-Hausdorff formula we obtain the following relationship 

W{a) = e"''-"^* = e-"''*e"*e-^. (2.1.4) 
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From canonical commutation relations f l2.1.ip we also find that for any /3, 7 G C 



el^^*h ={h- 13)6^^* and e^V = {b* + 7)e^^ (2.1.5) 
Denote the Weyl CCR algebra of observables generated by all the Weyl operators 

as 

l{{n) := {W{a) : a G C}. 
The linear map T : U{l-i) hl{'H) is called quasi-free, if it has the form: 

T{W{f)) = <il{f)W{T{f)) fen, (2.1.6) 

for some linear map F : H — ?■ H and complex function : H — C. 

2.1.3 States 

The state u on the algebra of observables U maps each observable A into its expec- 
tation value, which is in general a complex number u}{A) with the properties 

1) normalization: uj{1) = 1 

2) linearity: for each pair A, B of observables and each pair A, yU of complex numbers 
one has (jj{\A + ^B) = Xco^A) + fiuj{B) 

3) positivity: for each observable A, u{A*A) > 0. 

Let u be any state on the Weyl algebra hl{'H). This state is well-defined if for all 
/ G H all the expectation values uj{W{f)) are known, in other words it is completely 
defined by its characteristic functional: 



H3f^u{WU)). 
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We consider states such that their characteristic functional could be written in terms 
of correlation functions 

oo 

uj{W{f)) = u;{e^'^f^) = J2'-i.iKm 

n=0 

oo .„ 

= ^MY.h'^wr)t}, (2.1.7) 

n=l 

where the truncated correlation functions uj{...)t are defined recursively 

uj{b{h)...b{U)) = J2^{b{h)...),...oo{...b{U)),, 

here the sum is taken over all possible ordered partitions (1, ...), ...(..., n) of the set 
{1, n}. The expression uj{b{fi)...b{fn))t is called the truncated correlation function 
of order n. See [22] for more details. 

A state u of the boson algebra of observables is called a quasi-free state, if all its 
truncated correlation functions of orders n > 2 vanish. 

From (12.1.71) it follows that the general quasi-free state is determined by its one- 
and two-point correlation functions and therefore it is of the following form 

u{W{f)) = exp{za;(6(/)) - la;(6(/)6(/))i} 

= exp{zc.(6(/)) - liujibim - uj{b{f)r)} (2.1.8) 
= exp2r(/)-^s(/,/), (2.1.9) 



where we denoted r(/) = uj{b{f)) and s(/, /) = u;(6(/)^) — Lo{b{f)y. 

A state u is called regular if the function a t— )■ uj{W{af)) is continuous for all 
/ G "H and every a G C. Characteristic functionals of regular states on CCR(i3) are 
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characterized by Araki and Segal [H |2] in the following theorem. 

Theorem 1. A map 3 f ^ u){W{f)) G C zs the characteristic functional of a 
regular state oj on CCR(H) if and only if 

1. uiWifd)) = 1. 

2. The function a h-> uiW^af)) is continuous for all f ^Ti, a E C. 

3. For any integer n > 2, all /i, .../„ G "H and all zi, ...Zn G C one has 

n 

J2 ^(WU) - fk))e"'^'-^-f'''>/'z-z, > 0. 
j,k=i 

The Gibbs states or the equilibrium states are given by 

, Tr(e-'^'^*M) 

for any /3 G C. The state that is not in this form is called a non- equilibrium state. 
It can be shown that every Gibbs state of the form (12.1.101) is quasi-free (see [22] for 
details). 

2.2 Dynamics of quantum systems 
2.2.1 Reversible dynamics 

In quantum mechanics for a physical system described in terms of a Hilbert space 
"H a state G 7/ of the systems evolves in time under the action of a one-parameter 
strongly continuous group of unitary operators. 
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A family of unitary maps Ut G BiTi), t G M is called a strongly continuous one- 
parameter group of unitary operators if 

1) UtUs = Ut+sJoT s,tGM, 

2) \imt^o\\Ut(l)-(P\\=OJoi (Pen, 

3) ur' = U.f 

If {Ut}t<m. is a strongly continuous one-parameter group of unitary operators, then 
by Stone's theorem [] there exists a self-adjoint operator if, the Hamiltonian of the 
system, such that 

U, = e-*^, 

for t G M. 

The dynamics of the system can be described using two equivalent perspectives. 
Either a state changes in time while all obervables remain the same, of the dynamical 
group acts on the set of observables leaving the states unchanged. 

The second case is called "the Heisenberg picture". The time evolution of the 
observable A G B{'H) is described by 

Tt{A) = U*AUt = e^'^Ae"'". 



2.2.2 Lieb-Robinson bounds for reversible dynamics 

For Lieb-Robinson bounds the quantum systems are considered on the countable set 
of vertices F (called lattice) which is equipped with a metric d. A Hilbert space Tix 
is assigned to each vertex x G F. For any finite subset A C F the Hilbert space 



2.2. Dynamics of quantum systems 



19 



associated with it is the tensor product 



(2.2.1) 



The local algebra of observables over A is 



where B{l-Lx) denotes the algebra of bounded linear operators on l-Lx- 

If Ai C A2, then we may identify Ai^^ in a natural way with the subalgebra A\-^ 
11a2\Ai of Ak^, and simply write C Ap^.^. Then the algebra of local observables is 
defined as an inductive limit 



See [] for more detailed on the inductive limit. The C*-algebra of quasi-local observ- 
ables is the norm completion of A^^^. 

The support of the observable A G is the minimal set X C A for which 
A = A' ® 1a\x for some A' e Ax- 

We assume that there exists a non-increasing function F : [0, 00) — > (0, 00) such 
that: 

i) F is uniformly integrable over F, i.e., 




(2.2.2) 



Acr 




(2.2.3) 



and 
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ii) F satisfies 

Example 1. As an example one may take T = 7/ for some integer u > 1 with the 
metric d{x, y) = \x — y\ = Yl'j=i ~ Vjl- ^^^■^ case, the function F can be chosen 
as F{\x\) = (1 + for any e > 0. To show that the constant C is finite we use 

the triangle inequality and the symmetry of the taken supremum over x and y 



(1 + 


\x-y\Y+' 


(1 + 


\x - z\y+^{i + 


\y-z\y+- 



C = sup 

^ {l + \x - z\ + \y - z\ + iy+' 



< sup 



^ \1 + \x — z\ 1 + \y — z\ 



Now we use the inequality between geometric and algebraic mean: for any a > 1 and 
any a,b > 0, 



a + b\°' ^ a"' + b"' 



Continuing the calculations 



C< sup V2— 

VVl+|a:-2|/ \l + \y-z\J I 

= 2— '(sup V f-^)"' + sup v (-J—r-). 



Since two supremums are the same and each of them is achieved at any value we 



2.2. Dynamics of quantum systems 



21 



have 




1 



) 



1 + Ix — 2; 



< 00. 



Having a set F with a function F that satisfies f l2.2.3p and f l2.2.4p . we can define 
for any /i > the function 



which then also satisfies i) and ii) with ||-F)t|| < and < C. 

The Hamiltonian of the system is described by local Hamiltonians if'*"^ = {H^}, 
where is a not necessarily bounded self-adjoint operator over Hx, and bounded 
perturbations defined in terms of interaction $. The interaction $ is a map from 
the set of subsets of F to with the properties that for each finite set X C F, 
$(X) G Ax and = For any fi > 0, denote by i3^(F) the set of 

interactions for which 



F^{d)=e->^^F{d), 



(2.2.5) 



||$||^ = sup J2 



< 00. 



(2.2.6) 



For every finite subset A C F the Hamiltonian is of the form 



Since these operators are self-adjoint they generate a Heisenberg dynamics 
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for any A G ^a. 

The following theorem provides the Lieb-Robinson bounds. 

Theorem 2. Let X and Y be two disjoint subsets of A. Then for any pair of local 
observables A G Ax and B G Ay one has that 

\\[rt{A),B]\\ < ^^^^^(e^"*"''^''!*! -1) J] (2.2.7) 

xeXyeY 

The bound in the theorem could be rewritten using the properties of the function 
F as follows 

||[r/^(A),5]|| < M^iM||F||min(|X|,|F|)e-'^('^(^'^)-^-l*l), (2.2.8) 

2||$|| C 

where Vf, = " '"^^ is the Lieb-Robinson velocity. 

This theorem is a special case of the general Lieb-Robinson bound fll.2.3p that 
we are going to consider here, so the proof of this theorem is included in the proof in 
section 13.31 as a particular case and it also can be found separately in [17] . 

2.2.3 Existence of the dynamics in the thermodynamic Umit 

We assume that the set F, which is still equipped with a metric d, as a countable set 
with infinite cardinality. One can take, for example, F = Z*^, for some u > 1. 

The thermodynamic limit is taken over an increasing exhausting sequence of finite 
subsets An, C F. 

Theorem 3. Let /i > and $ G i3^(F). The dynamics Tt corresponding to $ exists 
as a strongly continuous one-parameter group of automorphisms on Ar such that for 
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allteR 



lim Wrl'iA) -Tt{A)\\ = 



for all Ae Ar- 

The proof of this theorem can be found in [T7j or, as a particular case, in section 
13.41 where we prove the existence of the thermodynamic hmit for the irreversible 
dynamics. 

2.2.4 Irreversible dynamics 

In general, if we consider an open system taking into the account the interaction 
between the system and an environment we have to consider a non-Hamiltonian 
system. The dynamical maps 7t of such system form a one-parameter completely 
positive dynamical semigroup on an algebra of observables U, which is described by 
the properties 

1) 7t is completely positive, 

2) Itil) = /, 

3) islt = 7t+s, 

4) 7t — > 11 when t — > strongly, i.e. for any A ^U, lim^^oo ||7t(^) — ^|| = 0. 

A map 7 is called completely positive if for any n > 1, the linear maps '-f^idM„ where 
Mn are n x n complex matrices, are positive. 

A strongly continuous group of automorphisms is a special case of a dynamical 
semigroup for 7t(A) = Tt{A). 

The dynamics 7t of the open system is the solution to the Markovian master 
equation 



d_ 
di 



MA) = C{j,{A)), 



(2.2.9) 
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for all t G M, with 70 (^) = A, where the generator C acts on the space of observables. 
It was shown in [13] that the generator C is in the Lindblad-Kossakowski form 



where La are such that J2a ^l^a is a bounded operator and here {A, B} = AB+BA is 
the anti-commutator. The Hamiltonian H have bounded interaction terms, but may 
have unbounded on-site terms as in section [2.2.2[ making the generator C unbounded. 

The proof that the time evolution of an open system is Markovian in the weak 
coupling limit can be found in [6]. It was shown there that the limit of the unitary 
dynamics that depends on a coupling constant may not be unitary, but it satisfies 
the Markovian master equation ( I2.2.9P when the coupling constant converges to zero 
for a rescaled time variable. 

In this dissertation we are going to discuss the dynamics generated by the observ- 
ables of the type (12.2. lOp . where H and each observable La are time-dependent. We 
will show that the evolution with this type of generator is well-defined as a norm- 
continuous cocycle of unit preserving completely positive maps. 

2.2.5 Quasi-free dynamics 

The unity-preserving quasi-free semigroup {Tj}j>o on the CCR algebra UiTi) is de- 
fined in a similar way: 




(2.2.10) 



a 



T,{W{f)) := %{f) W{T,{f)) , fen, 



(2.2.11) 
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where ^t=o(/) = 1, Tt=o = 1 and = 0) = 1, Ttif = 0) = . Note that the 

semigroup property of {Tt}t>o and (12.2. lip imply 

Ts+t{W{f))=T,{Tt{W{m (2.2.12) 
= v[/,(r,(/))vI/,(/) W{T,{T,{m = ^,^,{f)W{T,_,,{f)) . 

Then hnear independence of the Weyl operators yields 

r,+,(/) = T,{Tt{f)) and ^s+t{f) = ^s{Tt{f))%{f) • 

Hence, {rt}(>o is in turn a semigroup on Ti. 

Note that by definitions (12.2. lip and (I2.1.8P the quasi-free dynamics maps the 
quasi-free states into the states: 

Ur,s{Tt{Wim = ^tif) cOr,s{W{Ttim = ^*(/) ^n,sAW{f)) , (2.2.13) 



where n(/) := r{Tt{f)) and St{fJ) := s{Tt{f),Tt{f)). In general the states (12X131) 
are not quasi-free. 
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Chapter 3 

Lieb-Robinson bounds and the 
existence of the thermodynamic 
limit for a class of irreversible 
quantum dynamics 

3.1 Set up 

The set up here is the similar to the set up for the Lieb-Robinson bounds for uni- 
tary dynamics as in section I2.2.2I The system is considered on the countable set F 
equipped with a metric d. We put the same restriction on the lattice: there is a 
non- increasing function F : [0, oo] — )■ (0, oo) that satisfies fl2.2.3p and fl2.2.4p . 

The Hilbert space of states 'H\ of any finite subsystem A C F is defined by (12.2.11) 
as a tensor product of Hilbert spaces l-Lx of every point x in A. The C*-algebra of 
quasi-local observables is the norm completion of the algebra of local observables 
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A'^'^ defined in (12.2.21) as a space of bounded linear operators on Ha. 

Tlie generator of the dynamics is now consists of two parts: the Hamiltonian 
interactions and the dissipative terms. We allow both of them to be time-dependent. 

As in the case of unitary dynamics, the Hamiltonian part is described by an 
interaction ■), which is now time- dependent. For every time t G M the interaction 
■) is a map from the set of finite subsets of F to such that for every finite 
Z CT 

1. <^{t,Z)eAz, 



2. = <l>{t,Z). 

The dissipative part is described by terms of Lindblad form. For every finite Z cT 
these terms are defined by a set of operators La(t, Z) G Az, where a = 1, N{Z). 
It is possible for N{Z) = oo with an additional assumption of the convergence of the 
sum. 

Then, for any finite set A C F and time t G M the generator C\ G B{Aa,Aa) is 
defined as follows: for all A G Aa, 

£A{t){A) = J2 ^^zmA), where (3.1.1) 

ZcA 

^zmA)=^mt,Z),A] (3.1.2) 

N{Z) 

+ J2 [Llit, Z)ALa{t, Z) - -{La{t, zyiait, Z), A} j , 

here {A, B} = AB + BA is the anticommutator of A and B. The operator z{t) can 
be viewed as a bounded linear transformation on Aa for any A D Z, which can be 
written in the form z{t) ® id_4^ . For every Z C A the norm of these maps can be 
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bounded independently of the choice of A D Z as follows: 



N{Z) 



^zm<m{t,z)\\ + 2Y,\\La(t,z)\\\ 



a=l 



If N{Z) = oo we can insure that the sums Xlaii converge guaranteeing 

the uniform boundedness of the maps ^z(i)- However, it is more general to assume 
that the maps "^zit) defined on Az are completely bounded. A map \1/ G B{Az) is 
called completely hounded if for all n > 1, the linear maps \1/ ® idM„! where M„ = 
i3(C") are n x n complex matrices, defined on Az ® M„ are bounded with uniformly 
bounded norm. The cb-norm is then defined by 



By this definition the cb-norm of G B{Az,Az) is independent of any A such 

that Z C A C r. 

The main assumption that we make is the following. 

Assumption 1. Given {T,d) and F as described at the beginning of this section, the 
following hypotheses hold: 

1. For all finite A C F, £a(^) is norm- continuous in t (with the uniform operator 
norm on B{Aa, Aa)), and hence uniformly continuous on compact intervals. 

2. There exists fi > such that for every t G M 



^llcb = sup 11^ (8)idA/„|| < oo. 



n>l 




(3.1.3) 



where 



cb denotes the cb-norm of completely bounded maps. 
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The second part of the assumption means that the interaction weakens exponen- 
tially as the diameter of the support grows. This condition is similar to f l2.2.6p in 
unitary dynamics case. 

From definitions note that we have 

wum < w^zim <Y.Y. w^^mu < iiv&ii,,iaiiifii. 

ZCA x,yGAZ3x,y 

Define 

Mt=\\^\\tJA\\\F\\. (3.1.4) 

From the definition fl3.1.3p it is clear that Mg < Mt for s < t. 

To define the dynamics of the system fix T > and, for all A G Aa, let A(t),t G 
[0,T] be a solution of the initial value problem 

l-A(t) = £A(t)A(t), A(0)=A (3.1.5) 
at 

Since ||£a(^)|| < Mt < oo, this solution exists and is unique by the standard existence 
and uniqueness results for ordinary differential equations. 

For < s < t < T, define the family of maps {7t^s}o<s<t C B{Aa,Aa) by 

where A{t) is the unique solution of (13.1.51) for t G [s, T] with initial condition A{s) = 
A. 

Then, the cocycle property, '~ft,s{A{s)) = A{t), follows from the uniqueness of the 
solution of (I3.1.5p . 

Recall that a linear map ^ : A ^ B, where A and B are C*-algebras is called 
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completely positive if the maps 7 ® id : ^ ® M„ — )■ i3 M„ are positive for all n > 1. 

Here M„ stands for the nxn matrices with complex entries, and positive means that 

positive elements (i.e., elements of the form A* A) are mapped into positive elements. 

3.2 Existence of a dynamics of the finite system as 
a semigroup of completely positive unit pre- 
serving maps 

In this section we show that the dynamics defined in the previous section exists as a 
norm-continuous cocycle of unit preserving completely positive maps. This extends 
the well-known result for time-independent generators of Lindblad form [?] to the 
time-dependent case. The theorem is formulated for more general generators of the 
Markovian dynamics, which includes the generators of the type fl3.1.ip . 

Theorem 4. Let A be a C* -algebra, T > 0, and for t G [0,T], let C{t) be a norm- 
continuous family of bounded linear operators on A. If 
(^)C{t){l) = 0■, 

(ii) for all AeA, C{t){A*) = C{t){A)*; 

(m) for all A e A, C{t){A* A) - C{t){A*)A - A*C{t){A) > 0; 

then the maps ^t,s, ^ s < t < T , defined by equation l^3.1.5\} . are a norm- continuous 
cocycle of unit preserving completely positive maps. 

It is straightforward to check that the C\{t) defined in (I3.1.ip satisfy properties 
(i) and (ii). Property (iii), which is called complete dissipativity, follows immediately 
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from the observation 

N{Z) 

U{t){A*A) - U{t){A*)A- A*U{t){A) = J2^A,L,{t,Z)]*[A,L^{t,Z)] > 0. 

ZcA a=l 

Therefore, using this result, we conclude that, under Assumption [H for all finite 
A C r, the maps 7^^, < s < t, form a norm-continuous cocycle of completely 
positive and unit preserving maps. 

Before we begin the proof we make the simplification of some notations. 

Let C{t), t > 0, denote a family of operators on a C*-algebra A satisfying the 
assumptions of Theorem H] and for < s < t consider the maps A 3 A '-ft,s{A) 
defined by the solutions of (13.1.51) with initial condition A aX t = s. Without loss of 
generality we can assume s = in the proof of the theorem because, if we denote 
C{t) = C{t + s), then 'jt,s = ^t-s,o, where 74,0 is the maps determined by the generators 

The maps 74,5 satisfy the equation 

Jt,s = id + [ C{T)^r,sdT. (3.2.1) 
J s 

In our proof of the complete positivity of 74,0 we will use an expression for 'jt,o as the 
limit of an Euler product, i.e. approximations T„(t) defined by 

Ut) = f\(id+-Ci-)] . (3.2.2) 

k=n ^ ^ 

The product is taken in the order so that the factor with /c = 1 is on the right. 
Lemma 1. Let C{t), t > 0, denote a family of operators on a C* -algebra A satisfying 
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the assumptions of Theorem^ Then, uniformly for all t G [0,T], 



lim \\Tn(t) - 7f,o|| = 



where Tn(t) is defined by h3.2.2\) . 



Proof. From the cocycle property established in Section 13.11 we have 



Itfi = \{lt 

k=n 



k_ ^ , 



Now, consider the difference 



w - = n (id - n 



k=n 

n 



k=n 



n(id+^/:(t^)) 



k=n 



n ' n 



To estimate the norm of this difference we look at each factor separately. 

Using the boundedness of C{t) and the fact that M^, defined in (13.1. 4p . is increas- 
ing in t, the norm of the first factor is bounded from above by 



n(id+^£(t^))ii < n(i + < (1 + -M,r. 

k=n k=n 



To bound the second factor notice that from fl3.2.ip we obtain 



7mII < 1+ / \mr)\\\\jr,s\\dT. 



Then by Gronwall inequality [TT| Theorem 2.25] we have the following bound for 
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the norm of the 74 s- 



\\lt,s\\ <ei'*II^MII'^^ <e^^*(*-^). 



Using again (13.2.11) we can rewrite the second factor as follows: 



id+-£(t^)^ -7,. = -C{t^-^- 



Lit 



^s)-f^,j^ ds 



n 
n 

3-1 



n 



(7, 



id) 



C{t^—l)-C{s)]-C{s){x,^-\d) 



ds 



C{t^^—^)-C{s)]ds- [\c{s) I C{T)-f^^i^dTds. 



n 



Therefore, the second factor is bounded from above by 



id+-£(t^' 

n n ' 



n n 



n 



eir-t'^W^dTds 



(s~t 

j-1 n 



'-)ds 



n V 2n 



where e„— J-Oast/n— j-O due to the uniform continuity of C{t) on the interval [0,t]. 
The third factor can be estimated in a similar way: 



ii7,^,oii= n HKt^\\= n 

k=j—l k=j—l 



n 



fc=i-l ^ 

< (1 + -M,)'^. 
n 
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Therefore, combining all these estimates we obtain 



\Tn{t) - 7*,o|| < n{l + -M,r- ( 6„ + M^^^^V ) (1 + -M,y 

n n \ In J n 

2n / 



This bound vanishes as n — oo. □ 

To prove Theorem H] we use the Euler-type approximation established in Lemma 
[H We show that the action of T„(t) on a positive operator is bounded from below 
and this bound vanishes as n goes to oo. 

Proof of Theorem [H First, we look at the each term in the Euler approxi- 
mation Tn{t) separately. For any t and s the complete dissipativity property (iii) of 
£(s), assumed in the statement of the theorem, implies 

< {id+tC{s)){A*){id+tC{s)){A) = {A* + tC{s){A*)){A + tC{s){A)) 
= A* A + tA*C{s){A) + tC{s){A*)A + t^C{s){A*)C{s){A) 
< A* A + tC{s){A*A) + t^C{s){A*)C{s){A). 

Since {C{s){A)y{C{s){A)) < \\C{s)f\\A\\, one gets 

< {id+tC{s)){A*A)+t^\\C{s)f\\A\\^ (3.2.3) 
< {id +tC{s)){A* A) +t^M^\\A\\^. (3.2.4) 

Let us apply the above inequality to the operator!?, where B*B := — A*A. 
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Note that \\B*B\\ < , so \\B\\ < \\A\\. 

< {id +tC{s)) {\\Af - A* A) +fM^\\Af (3.2.5) 

= \\Af - {id+tC{s)){A*A)+t^M^\\Af (3.2.6) 

From the (l3X3|) and (13:231) we obtain 

-t^M^WAf < {id +tC{s)){A* A) < {l+t^M^)\\Af (3.2.7) 

and therefore: 

-{l+t^M^)\\Af < {id +tC{s)) {A* A) < {l + t^M^)\\A\\\ 

So we get 

\\{id+tC{s)){A*A)\\ < {l + t'M^)\\Af . (3.2.8) 

Now, in order to bound the approximation T„(t) we first derive the following 
auxiliary estimate. For any fixed n > 1 we have: 

1 ^ n-l 

l[{id+sC{ks)){A*A) > -s'\\ArMl{l + -)"-l5^D(.)^ (3.2.9) 

J. J. n — 1 ^ — ^ 

k=n k=0 

where the value of s is chosen to be such that 



D{s) := 1 + s'Ml < (1 + ^)"-V(l + -^r-', (3.2.10) 

n — 1 n — 2 



with the convention that (1 H — ""^ = 1, for n = 1. 

We prove this claim by induction. The statement holds for n = 1 by f l3.2.5p . Now, 
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assume that fl3.2.9|] holds for n — 1. Then 

1 ^ n-2 

n {id+sc{ks))iA*A)+s'\\ArMi_,^,{i + ^r-'Y.^(')' ^ 

k=n-l k=0 

Since the left-hand side is a positive operator, we can write it as B*B. Then, 
1 

Y[{id+sC{ks)){A*A) = {id +sC{ns)){B*B) 

k=n 

n-2 

- s'\\ArMi_,),{i + ;^r'E^(^)' 

k=0 

^ n-2 

> -s'Ml\\B*B\\ - s'WAfMlil + E ^(^)' • 

k=0 

Here, we used fl3.2.5p and the fact that Mt is monotone increasing. This gives the 
following upper bound for 

n-2 

k 



\\B*B\\< n m+sC{ks)){A*A)\\+s^ArMl_,^,{l + -^r^^ 

k=n-l k=0 
1 ^ n-2 

< n i^+^"Mi)\\Ar+s'\\ArMi{i+^r-'j2D{s 



1 , n-2 

n-2' 

k=n-l k=0 

1 ^ n-2 

k=n-l k=0 



^ n-2 



k=0 
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Therefore we obtain 



Yl {l + sC{ks)){A*A) 

k=n—l 

^ n-2 

> -s'MlWAfDis}^-' - s'Mlis'Ml + 1)(1 + ^ 

n — 2 ^-^ 

k=0 

n — 1 n — I ^ — ' 

fc=0 

^ n—l 

>-s'Mi{i+ — -r~'\\Arj2D{s)\ 

n — 1 ^ — ^ 



fc=0 



where to pass to the second inequality we use our assumption on s (13.2. lOp . This 
completes the proof of the bound fl3.2.9l) . 

To finish the proof of the theorem we use Lemma [1] to approximate the propagator 
and put s = ^ in the bound fl3.2.9p . which yields 



n-1 



+ -C{-)){A*A) > --\\ArM^{l + (3-2-11) 
■^■^ n n n — l n 

k=n k=0 

Since I^l^)" = (1 + ^M^Y 1 as n cx), we get the estimate D{^^Y < 2 for 
1 < k < n. The right hand side of (13.2. lip is bounded from below by — ^||/l|pe M^2n, 
which vanishes in the limit n oo. 

To show the complete positivity of 74^0 note that any generator C\{t) satisfying 
the assumptions of the theorem can be considered as the generator for a dynamics on 
A B{C"'), for any n > 1, which satisfies the same properties, and which generates 
7(^cj®id acting on A®B{C"-). By the arguments given above, these maps are positive 
for all n. Hence, the 74^^ are completely positive. □ 
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3.3 Lieb-Robinson bounds for a class of irreversible 
dynamics 

For X C A, let Bx denote the subspace of B{Ax) consisting of all completely bounded 
linear maps that vanish on 1. 

It is important for us that all operators of the form 

^ 1 
JCxiB) := z[A, B] + Y^iKBLa - -{LIL,, B}) , 

a=l 

where A, La G Ax, belong to Bx, with 

N 

\\}CxU<2\\A\\+2j2\\Laf- 

a=l 

In particular, operators of the form [A, ■] appearing in the standard Lieb-Robinson 
bound fl2X7D special case of this general form. 

We can regard ICx as a linear transformation on Az, for all Z such that X G Z, 
by tensoring it with id^^^^^ , and all these maps will be bounded with norm less then 

ll^xllcb- 

Theorem 5. Suppose Assumption{^ holds. Then the maps satisfy the following 
bound. For X,Y G A, and any operators /C G Bx and B G Ay we have that 

^ xeXcAyGYcA 

The proof of the Lieb-Robinson bounds for 7^^ is based on a generalization of the 
strategy [15] for reversible dynamics, and on [20] for irreversible dynamics with time- 
independent generators. This allows us to cover the case of irreversible dynamics 
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with time-dependent generators. 

Proof of Theorem [5} Consider the function / : [s, oo) ^ defined by 

fit) = /C7t(5), 

where /C G Bx and B G Ay, as in the statement of the theorem. For X C A, let 
X"^ = A \ X and define Cx'^ and Cx by 

Cx^it) = Yl ^^(^) 

Cxit) = Cxit)- Cx^it). 



Clearly, [/C, Cx<^{t)] = 0. Using this property, we easily derive the following expression 
for the derivative of /: 

f\t) = Jcc{tht(B) 

= Cx4t)}C^t{B) + ICCxithtiB) 
= Cx4t)f{t) + }CCx{thtiB) , 

Let 'j^^ be the cocycle generated by Cx<^{t)- Then, using the expression for f'{t) we 
find 

f{t) = l^:m + [\^:iCCx{rhUB)dr . 



Since 7^/ is norm-contracting and ||/C||cb is an upper bound for the ||/C|| regarded as 
an operator on ^a, for all A, we obtain 

wmw < wmw + mu f \\cx{r)^tm\dr. (3.3.1) 
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Let us define the quantity 

CB[X,t) := sup — — . 

r^Bx II / II cb 

Note that we use the norm ||T||cb, because, as mentioned before and in contrast to 
the usual operator norm, it is independent of A. Then, we have the following obvious 
estimate: 

CBiX,s) < \\B\\6YiX), 

where 5y(X) = OifXny = and 5y(X) = 1 otherwise. From the definition of 
the space Bx we get that T(-B) = 0, when T G Bx, since B has a support in Y and 

Fnx = 0. 

Therefore (13.3. ip implies that 

CB{X,t)<CB{X,s)+ f \\Cz{r)\\CBiZ,r)dr. 

Iterating this inequality we find the estimate: 

°° (j- \n 

CB(X,t)<||i?||V^^^a„, 

n=0 

where: 

for n > 1 and = 1, (recall that is a constant, that appears in a definition of 
F^). The following bound immediately follows from this estimate: 
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If we take an exponentially decaying function F, defined in f l2.2.5p as Ff^{d) = 
e~^'^F[d) we can rewrite the above bound as follows 

11^7^.(5)11 < M£iM||F|| niin(|X|, |y|)e-^W^'^)-^^(*-)). (3.3.2) 

So the Lieb- Robinson velocity of the propagation for every t G R is 

vt,t, ■= • 

□ 

Note that the Lieb-Robinson bound (13.3.21) depends only on the smallest of the 
supports of the two observables. Therefore one can get a non-trivial bound even when 
one of the observables has finite support but the support of the other is of infinite 
size (e.g., say half the system). 

The bound in the LIbe-Robinson bound is uniform in A. This is important for 
the proof of existence of the thermodynamic limit of the dynamics, which is one of 
the main applications of Lieb-Robinson bounds, which will be presented in the next 
section. 



3.4 Existence of the thermodynamic limit 

The setup for the analysis of the thermodynamic limit can be formulated as follows. 
Let r be an infinite set such as, e.g., the hypercubic lattice Z'". We prove the existence 
of the thermodynamic limit for an increasing and exhausting sequence of finite subsets 
A„ C F, n > 1, by showing that for each A G Ax, {lt[s{^))n>i is a Cauchy sequence 
in the norm of Ar- To this end we have to suppose that Assumption [1] (2) holds 
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uniformly for all A„, i.e., we can replace A in f l3.1.3p by F. 

Theorem 6. Suppose that Assumption^ holds and, in addition, that Ii3.1.3\) holds for 
A = r. Then, there exists a strongly continuous cocycle of unit-preserving completely 
positive maps 'jf^ on Ar such that for all < s < t, and any increasing exhausting 
sequence of finite subsets A„ C T, we have 

hm ||7t"(^)-7M(^)ll=0, (3.4.1) 

for all A e Ar- 

The proof of existence of the thermodynamic limit mimics the method given in 
the paper [15]. 

Proof of Theorem [6j Denote £„ = Ca„ and 7^^; = 7!"^ Let n > m, then 
Am C A.„ since we have the exhausting sequence of subsets in T. We will prove that 
for every observable A G Ax the sequence {■y^g{A))n>i is a Cauchy sequence. In 
order to do that for any local observable A G Ax we consider the function 

fit) :=7ff(^)-7l?(^)- 
Calculating the derivative, we obtain 

f{t) = C^^i:^{A)~Cmli:\A) 

= £„(t)(7ff(A) - 7;? (A)) + (£„(t) - Cmmi-\A) 

= Cnit)f{t) + (c^it) - cut)hi::\A). 
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The solution to this differential equation is 

J s 

Since 74^'^) is norm-contracting, from this formula we get the estimate: 



ll/WII< / WiUr) - CUr)H::\A)\\dr 

J s 

^ E E /Vz(r)(7(?(A))|Mr. 

Using the Lieb- Robinson bound and the exponential decay condition ( ]3.1.3p . which 
we assumed holds uniformly in A, we find that 

ft 



zeA„\A„i x6X 

ft 



< Pillion /" e'^^'-'''('-^)rfr|X I sup F{d{x,z)). 

ZGAn\Am 

Since F is uniformly integrable fl2.2.3p the tail sum above goes to zero as n, m — )■ cxd. 
Thus 

ll(7ff-7l?)(^)ll^0, as n,m^ 00. 

Therefore the sequence {'^t"s {A)}'^=o is Cauchy and hence convergent. Denote the 
limit, and its extension to ^r, as 7^^. 
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To show that is strongly continuous we notice that for < s < t,r < T, and 
any A G Ay'^, we have 

ll7L(^) -7^.(^)11 < ll7L(^) -7^(^)11 + Il7;?(^) -7^J(^)II + ll7^J(^) -7.^.(^)11, 

for any n G N such that A G Aa„. 

The strong continuity then follows from the strong convergence of 7^"'' to ^yf^, 
uniformly in s < t G [0,T], and the strong continuity of 7f"'* in t. The continuity of 
the extension of 7]"^ to all of A G At follows by the standard density argument. The 
argument for continuity in the second variable, s, is similar. □ 
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Chapter 4 

Non-equilibrium state of a leaking 
photon cavity pumped by a 
random atomic beam 

4.1 Set up 

4.1.1 Description of the model 

Our model consists of a beam of two-level atoms that passe one-by-one a photon 
cavity C. Atoms in a beam are randomly excited. During the passage time r the 
corresponding single atom is able to interact with the cavity. For simplicity we 
consider a, so called, tuned case, when the cavity size is equal to the interatomic 
distance, so there is always one atom in the cavity. 

The cavity is a one-mode resonator described by quantum harmonic oscillator 
with the Hamiltonian He = €b*b ® 1 in the Hilbert space He, where b* and b stand 
for boson (photon) creation and annihilation operators with canonical commutation 
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relations (CCR): [b,b*] = 1, [b,b] = [b*,b*] = 0. 

The beam of two-level, {E,0}, atoms is described as a chain = X]n>i -^^n 
of individual atoms with Hamiltonian if_4^ = 1 ^ Erjn in the Hilbert space = 
®n>i'HAn- Here for any n > 1, T-Lau = the individual atomic operator 

rjn '■= (cr^ + l)/2, where is the third Pauli matrix. The eigenvectors f]n'^n ~ '^t, 
and Tjnip^ = 0, are interpreted as the excited and the ground states of the atom, 
respectively. 

The initial state of the system is the product state of the cavity and the states of 
each individual atom: 

Ps ■■= PC ® 6d Pk ■ (4.1.1) 



k>l 

Here pc is the initial state of the cavity, which we assume to be normal, i.e., given by 
a density matrix pc G CiiTic), where (Ti denotes the space of the trace-class operators 
on Tic, and ■= ®k>iPk is the state of the atomic beam. The states {pA:}fc>i on 
the algebras {Ci('H^j.)}fe>i are diagonal and identical (homogeneous beam), hence of 
the form 



pk 



p 



Here, the parameter p := Tr-^^^ (r/fc Pk) denotes the probability with which each atom 
is in its excited state 

In our model the repeated cavity-atom interaction is time-depended and piece- 
wise constant, that has the form: 

Knit) X[(n— l)T,nT) 

(t) A (6* + 6) ® r/„, . (4.1.2) 



Here Xiip^) is characteristic function of the set /. 
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The Hamiltonian of the entire system acts on the space T-is '■= T^c ® 'Ha, and is 
given by the sum of the Hamiltonian of the cavity, of the atoms, and the interaction 
between them 

H{t) = Hc + J2{Ha,. + K^it)) (4.1.3) 

n>l 

X[{n— l)T,nT) 

{t){\{h* + h) (g)r]n) . 



n>l n>l 



Notice that for the time t G [{n — l)r, nr), only the n-th atom interacts with the 
cavity and the Hamiltonian in time- independent. 

4.1.2 Hamiltonian dynamics of perfect cavity 



Let t E [{n — l)r, nr). Then the Hamiltonian (14.1.31) for the n-th atom in the cavity 
takes the form 

Hr,:=eb*b0l + l®Er]^ + \{b* + b)®r]r, . (4.1.4) 

Even though the atomic beam is infinite, there is only a finite number of atoms 
that have interacted with the cavity at any finite time t. Assuming that the initial 
state of the system is normal, our system can be described by normal states ujs{-) : = 
Tr( ■ ps), which are defined by trace-class density matrices ps G CiiTic ® 'Ha)- 

Then the partial traces over T-La and over Tic: 

ujc{-) := Tr^^(- ps) and w^(-) := Tr^^(- ps) , (4.1.5) 



define respectively the cavity and the beam states. 

We often refer to density matrices as states, if this wording will not produce any 
confusion. 
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We suppose that initially our system is in a product state: 

ujs{-) := (ujc ® ujA)i-) = Tr( ■ pc ® Pa) , ps = pc ® pA , (4.1.6) 

where pc E £i('Hc) and G £1(7/^). 

For any state ps on the algebra of observables 2l('Hc ® "H^i) the Hamiltonian 
dynamics of the system is defined by (14.1.31) . or by the quantum time-dependent 
Liouvillian generator: 

Lit)ips):=-t[Hit),ps] . (4.1.7) 

Then the state ps{t) of the total system at the time t is a solution of the Cauchy 
problem corresponding to Liouville differential equation 

j^psit) = mipsit)) , ps{t = 0) := PC® PA- (4.1.8) 

We denote by ujg{-) := Tr( ■ psif)) the system time evolution due to (I4.1.8P for the 
initial product state ujs{-) (14.1.61) . 

Notice that in general the Hamiltonian evolution (I4.1.8P with time- dependent 
generator is a family of automorphisms {Ff <,}o<s<t with the cocyde property: 

Ps{t) = Tt,sPs{s) . (4.1.9) 

For our model with tuned repeated interactions the solution of (I4.1.8P and the form 
of the evolution operator (14.1.90 are considerably simplified. Indeed, the generator of 
the dynamics of our system (14.1.30 is time-independent for each interval [(n — l)r, nr). 
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Therefore, by virtue of fl4.1.4p and f l4.1.7p . the LiouviUian generators 
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Ln := L{t) , te[{n- l)r,nr) , n > 1 , (4.1.10) 

are time-independent and commuting. Note that any moment t > has the repre- 
sentation 

t ■=n{t)T + v{t) , n{t) = [t/r] and z/(t) G [0, r) , (4.1.11) 



where [x] denotes the integer part of the number a; G M. Then solution of (I4.1.8P for 
t G [{n — l)r, nr) gets the form: 

Ps{t) = Tt,s=o{pc ® Pa) = e^^"e^^"- ... e^^^e^^HPc ® Pa) ■ (4-1.12) 

In the next section we exploit a specific structure of the Hamihonian dynamics 
fl4.1.12p and a special form of interaction fl4.1.2p to work out an effective evolution 
of the perfect cavity Hq. Our results concern first of all the evolution of the photon 
number expectation in the time-dependent cavity state 

pcit) := Tin^psit) . (4.1.13) 

Note that by fl4.1.3p . fl4.1.4p one gets that in our model the atomic states pk do 
not evolve: 

[1® Pk,H{t)]=0 , k>l. (4.1.14) 
The form of the initial state (14. 1.1 p together with (I4.1.12p . (l4.1.13p imply a discrete 
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time evolution for the cavity state given by the following recursive formula: 

n 

pP ■=pc{t = nr) = Tiu^psinr) = Tr^^ [e^^"...e^^^e^^Hpc ® 0Pfc)] (4.1.15) 

A;=l 

n-1 

=Tr«^Je-^"{Tr«^_^...Tr«^^e-^"-...e-^^e-^nPc ® (g)pfc)} ® Pn] 

k=l 

For any density matrix p G (ti('Hc) corresponding to normal state on the operator 
algebra 2l('Hc) we define the mapping C : p ^ '^{p), by 

Cip) := Tr^^Je-^"(p® p„,)) = Tr«^Je"^-^"(p ® p„)e*-^"] . (4.1.16) 

Here the last equality is due to fl4.1.7p and (14.1.101) . 

Note that the mapping (14.1.161) does not depend on n > 1, since the atomic states 
{Pn}n>i are homogeneous. Then the cavity state at t = nr is defined by the n-th 
power of £: 

pP = C{pt'^) = C-{pc) . (4.1.17) 

Therefore, by (]4.1.12p . ( 14.1.15^ and (14.1.17P one obtains that for any time t = nr + 1^, 
where u G [0, r), the cavity state is 

pc{t) = Tr^,^^^ [e'^^"+H^"(Pc) ® Pn+i)] • (4.1.18) 

We are interested in the expectation value N{t) of the photon-number operator 
:= b*b in the cavity at the time t: 



N{t) ■.= oo'c{N) = TTnc{b*bpc{t)) 



(4.1.19) 
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For t = nr the state of the cavity can be expressed using fl4.1.17p . which gives 



51 



Ninr)=TTnAb*bC^ipc)) • (4.1.20) 
4.1.3 Quantum dynamics of leaking cavity 

For a more general situation we take into account a possible leakage of the cavity, 
where photons may leave the cavity at some non-zero rate a > 0. We consider 
this case in the framework of Kossakowski-Lindblad extension of the Hamiltonian 
Dynamics to irreversible Quantum Dynamics with time- dependent generator 

L4t){ps) = -i[H{t),ps]+ab{ps)b* - ^{b*b,ps} , (4.1.21) 

for any ps E CiCHc ® H^). 

Similar to fl4.1.8p the evolution of the state is defined by a solution of the non- 
autonomous Cauchy problem corresponding to time-dependent generator ( I4.1.2ip 

j^psit) = L^{t){psit)) , psit = 0) := PC® PA- (4.1.22) 

In general for a time-dependent generator Lcr(t) the proof of existence of this solution 
is a non-trivial problem, as in chapter [3l In the present case of the tuned repeated in- 
teractions, when the Hamiltonian is time-independent for each interval [{n — l)r, nr), 
the generator fl4.1.2ip gets the form: 

LaAps) ■■= -i[Hn, ps] + (^b{ps)b* - ps} ,ps = Pc®Pa- (4.1.23) 

With restriction to the tuned case and with the help of fl4.1.23p the solution of 
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fl4.1.22p for t G [(n — l)r, nr) becomes of the form: 

ps{t) = t/t%o(Pc ® pa) = e'^^-"e"^-"-^ ... e"^-^e^^-i(pc ® Pa) • (4.1.24) 

Here {f/fs}o<s<t is a family of unit preserving completely positive maps with a cocycle 
property f l4.1.9p . 

The discrete evolution operator fl4.1.16p has to be modified for cr > as follows: 
Definition 1. For any state p on 2l('Hc) we define the mapping 

C^ip) := rr^^Je-^-(p®p„)) • (4.1.25) 

Let Pc '■= Pc{t = 0) denote the initial state of the cavity. Then similar to fl4.1.18p 
we obtain for pc{t) at the moment t = nr + z/, where u G [0, r), 

pc{t) = Tr^^^^^ [e'^^-+HC(Pc) ® Pn+i)] . (4.1.26) 

Now we define the functional w^(') Tr^p( . pcif))- We study the infinite-time 
limit for the cavity state ujc{-) '■= limt_j.oo i^M'). Here the limit means trace- norm 
convergence of the sequence (14.1.260 to density matrix pc,a '■ 

lim ||pc,.-pc,<x(t)||i = 0, (4.1.27) 

where the norm || . ||i denotes the trace-norm on the space of the trace-class operators 

We consider the functional 

uc{W{a)) = lim uj'c{W{a)) , (4.1.28) 

t—>-oo 
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generated by the Weyl operator on T-ic'- 
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W{a) = e""-"" , aeC . 

See Section 12.1.21 for more on Weyl operators. Notice that convergence (14.1.281) for 
the family of the Weyl operators guarantees the weak-* limit [2] of the states ^ci') 
when t — )■ CX3. 

4.2 Number of photons in perfect cavity 

First we consider the case of the perfect cavity, i.e. a = 0. Then for the discrete 
evolution operator one gets C„=o = C, see (14.1.161) and (14.1.251) . 

Our first result concerns the expectation of the photon- number operator N = b*b 
in the cavity (I4.1.19p . For t = nr this expectation value takes the form (14.1.201) . 

In the theorem below we suppose that the initial cavity state wc(-) is also gauge 
invariant, i.e. e*°^pce~*°'^ = Pc- 

Theorem 7. Let pc be a gauge invariant state. Then for t = nr the expectation 
value ^.1.20^ of the photon number in the cavity is 



2A^ 2A^ 
N{t) =iV(0)+rap(l-p) — (l-coser)+p^— (1-cosner) . (4.2.1) 

^2 ^2 



If for the initial state pc one takes in the theorem the Gibbs state for photons at 
the inverse temperature /3: 

= e-^^'*VTrwce-^''*' , (4.2.2) 
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then the number of photons (14.2.11) is 



1 2A^ 2A^ 

N(t) = -3 h -p) — (1 - coser) — (1 - cos ner) . (4.2.3) 



To prepare the proof of Theorem [7] we calculate first explicit expressions for C{p) 
acting on the space p G ^{T-ic) of the cavity states and for the n-th power of the dual 
operator C* that we apply to the number operator 6*6, see (I4.2.10p and Remark [H 

We split these calculations into two lemmas. The proof of the Theorem [7] is 
presented at the end of this section. 

The Hamiltonian (I4.1.4p can be written in the following form 

Hn =e {h* + ^?7„)(6 + ^?7„) + Erfn - 
=e b*b +{E- AVe)r/„ • 

Here b := b + Xf]n/^ and b* := b* + A?7„/e are new boson operators since by the CAR 
and CCR properties of ?7* and respectively b, b* one gets: [b, b*] = [b, b*] = 1 and 
[6, b] = [b*, b*] = for any index n. 

This is in fact a unitary shift transformation generated by 



Vn:=-{b* -b)(^r]n. (4.2.4) 



Then the transformed Hamiltonian is easily calculated to be given by 

Hn = e'^'-Hne-'"'" = e b*b + {E - ^)r/„. (4.2.5) 
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From the definition of Vn we obtain that 



To compute b = e*^" (]l ® 6)e we note that if 



then 



Therefore 



du 



e"' [h\h]e 



F{v) = F(0) - u. (4.2.6) 
Applying this formula to b we find that 



b = b(g)l~l(S) -r]n, 
e 



and similarly 



b* = b* -rin- 
e 



Note that the dynamics generated by Hn (or by leaves the atomic operator 
rjn invariant 

Lemma 2. For any state p on T-Lq 

C{p) = pe-^(b''b)/^e-'^''*'e^^'*-'^/'pe-^^''-'^/'e'^''*'e^^''-''^^ (4.2.7) 
+ (1 - p)e-'"^''* V"^*'' • 
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Proof. Using the shift transformation f l4.2.4p C can be written in the form 

£(p) = Tin^^ [e-'^"e-'^^"e'^-ip ® pn)e-'^"e'^^"e'^"]. (4.2.8) 

From definition of Vn and the fact that ?7„ commutes with p^^^ we have 

e'^'-ip ® Pn)e-^^" = e^(6*-'')Ape-"(^*-^)/^ ® r^^p, + p ® (1 - Vn)Pn. (4.2.9) 

Therefore, plugging this expression into fl4.2.8p we obtain 

+ Tin^^ [e-'^"e-*^^" (p ® (1 - r/„)p„)e^^^"e'^"]. 
Using the diagonal form (14.2.50 of the Hamiltonian, we have 

£(p) =TTnA„ ^^~X{b^^b)/e^-^reb*b^X{b'-b)/ep^^Xib'^b)/e^ireb'b^X{b'-b)/e ^ 

+ TiH^Je-'^''''''e^^''*-^^/'pe"^^^*~''^/'e'^'^*^ ® (1 - Vn)VnPn] 

+ Tln^^ ^^~Xib^~b)/e^-^reb*b^^^r.b'b^X(b^-b)/e ^ ^^^^ _ 

+ Tr^^Je"^-"* V^^'*' ® (1 - ^n)Pn]. 
Since (1 - r]n)Vn = 0, we get 



C{p) = Tr«^Je"^(''*-^)/^e-^-^''*^e^(^*-'')/>e-^(^*-'')/^e^-^''*V(''*-^)/^ ® r/„p„] 

-ireb'b „ATtb*b 



TrH.„ [e-^^^" V^^" ' ® (1 - ^n)Pn] 



^ g-A(b- -b)/eg-ire6*fegA(fe* -6)/e^g-A(b* ~b) / e ^ireb'b ^Xib* -b)/e 

+ (1 -p) e-^"^^*V"''*^ • 
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□ 

To calculate the expectation of the photon- number operator N = b*b at t = nr, 
we would need to find the action of the n-th power C"'{p) of the operator fl4.1.16p . 
But in fact it is easier to calculate this mean value using the n-th power of the adjoint 
operator £*, which for any bounded operator B G B{'Hc) and p G ^liTic), is defined 
by relation 

TrneiC*iB)p) := Tr^,(5 Cip)) . (4.2.10) 

Using the property ( I4.1.14p one can obtain explicit an expression for the operator C*. 
Indeed, by (14.2.101) and by fl4.1.1(i|] 

TthAB C{p)) =lTnc<^H^A{B ® 1) e-^^^"(p ®t){t® p„)e^"^"} (4.2.11) 
=Tr^,««^J(]l ® Pn)e~"''-{B ® ]l)e'-^"(p ® t)} 
=Tr^,^^^Je^-^"(i? ® p„) e-*-^" (p ® 1)} , 

where we used the commutator fl4.1.14p and the cyclicity of the trace Tr-^^^^^^^ . Then 
04.2.111) yields 

C{B) = Tr«^Je^-^"(5 ® pje-^^^"} , (4.2.12) 
which according to f l4.1.16p . is independent of n > 1. 

Remark 1. For density matrix p G Ci(Hc) such that Try^^{V{h,h*) C{p)) < oo for 
any polynomial V{b,b*) , one can extendthe definition lj^4-2.10\ ) of C* to this class of 
unbounded observables 2l('Hc). The advantage of using the dual operator C* is that 
its consecutive application does not increase the degree of the polynomials generated 
by operators b and b* . 



4.2. Number of photons in perfect cavity 58 
Now following the line of reasoning of Lemma |2] we deduce from (14.2. 12p that 

= pe-^(''*-^)/^e^"'^*V(^*-'')/'5e-^(''*-'')/'e-*"'''*V(''*-'')/' (4.2.13) 
+ (1 - p)e-'^'^'''Be'^'^''' . 



Lemma 3. For B = b*b and for the adjoint operator C* defined by ^.2.13\ ) we obtain: 



{C*TiP*b) = b*b + p ^[(1 - e""") b* + {l- e"™^") b] (4.2.14) 
+ np{l — p) ^-(1 — coser) + p —^(1 — cos ner) . 



Proof. From (I4.2.6P one gets 

^x{b*-b)/e^#^-\{b*-b)/e = 5# _ x/e, where b* = b or b. 

And therefore 

eA(b*-fe)A^*^g-A(fe*-b)A _ _ x/e){b - A/e). 
From the CCR properties of b and b* we find that their evolution is 

^iTeb*bjj* ^—ireb'b g*^^5* 

Therefore, making the shift to calculate the first term in (14.2.131) we get 

e-Hb*-b)/e(^^iTe^* _ A/e)(e-^^^6 - A/e)e^(''*-^)/^ (4.2.15) 
= {e'^'b* - (1 - e'^')\/e){e-'^'b - (1 - e-'^')\/e). 
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Hence, f l4.2.ia|) and fl4.2.15p imply 

C*{b*b) =pe-^(''*-^)/^e^"^^*V(^*-'')/^6*6e-^(''*-^)/^e-^"^''*V(''*-^)/^ (4.2.16) 
+ (1 - p)e-'^'^'%*be'^'^*^ 

=p{e''^b* - (1 - e'^")A/e)(e-""6 - (1 - e~''^)X/e) + (1 - p)b*b 

A A 2A^ 

=b*b + - e''^)b* +p-{l- e-''^)b + - coser) . 

If in ( 14.2. 13p we put B = b*, then one gets 



Now we are going to prove the n-th power formula (14.2. 14p by induction. Note 
that if we take n = 1 in this formula we get (14.2. 16p . 

Suppose that (14.2. 14p is true for n to check that it is also valid for + 1. 



C*(b*) = e''^b* -p{l-e''^)\/e , 



(4.2.17) 



and similarly, one obtains: 



C*{b) = e-''^b - p{l - e-''^)\/e , 



(4.2.18) 



for B = b. 



{c*Y^\b*b) =c*{{c*Y{b*b)) 




A 



Taking onto account (I4.2.16p . (I4.2.17P and (14.2. ISp we can express the action of 
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operator C* as follows 



A,. ...... A,. 2A_2 

,2 



=b*b + p-{l - e''^)b*+p-{l - e~''^)b + p—{l - coser) 



e e 
A 

e 

A A^ 

e 

2A^ , , 22A^ , , 22A^, 

+ — coser) — p -—-n[l — coser) + p — — cos ner)) 

^2 ^2 ^2 

--b*b + p-^{l - e("+^)*^")6* + p-^{l - e-("+^)"^)6 

2A^ / \ / \ 2 A^ / 

+ p^T-[l — coser)(n + 1) — p -^[2 — 2 coser — 2 cos ner 

^2 ^2 

+ 2 cos(n + l)er + 2n — 2n cos er — 2 — 2 cos ner). 



Simplifying the last expression we get 



{CY+\b*b) =b*b + p-{l - e("+i)'^^)6* + p-{l - e-("+^)^^^)6 

2A^ 2A^ 
+ p^-{l — cos er)(n + 1) — P^^ri^ + ~ coser) 

2A^ 

+ ^ cos(n, + l)er), 

e^ 



which proves (14.2.141) formula. □ 

Proof, (of Theorem [7j) To find the number of photons in the cavity at time t = nr 
we take the expectation in (14.2. 14p . Since the initial state is gauge invariant we get 



N{t) =Tr^,(6*6£"(p^)) = Tr^,((/:*)"(6*6)p^ 



2A^ 2A^ 
=A^(0) + np{l — p) -—- (1 — cos er) + p^— — (1 — cos ner) 

^2 



□ 
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Remark 2. Theorem implies that only flux of randomly exited atoms (i.e. < p < 1) 
is able to produce a pumping of the cavity by photons. Since the number of photons 
in cavity unboundedly increases with time, there is no regular limiting state. 



4.3 Number of photons in leaking cavity 

In the Hamiltonian case the number of photons in the cavity increases in time. To 
stabihze the energy in the cavity we consider the case when the generator of the dy- 
namics ( 14.1.2ip has a nonzero dissipative part, which describes the leaking of photons 
out of the cavity. So suppose that a > 0. 

The next theorem shows that the number of photons in the leaking cavity stabilizes 
in time for any a > 0. It is obviously different to the case a = 0, see Theorem [71 

Theorem 8. For any a > and arbitrary gauge-invariant initial cavity state pc such 
that 

co'c{b*b) \t=o= Trnc{b*b pc) < oo , (4.3.1) 
one gets for the limit of the cavity photon-number mean value 



iV~ = cuc{b*b) := lim 4(6*6) = (4.3.2) 

t—^oo 

^ -{l + e-'^"(l -2p) -2e-'^"/2(l -p)coser} . 



4e2 + a2 1 - e-'^^ 
Lemma 4. For any state p on algebra QilHc) one has: 

C,(p) = pe-A{^--^)Ae-^?(e^(''*-^)/>e-^(^*-'')/^)e^(''*-^)/^ + (1 - p)e^^f (p), 



4.3. Number of photons in leaking cavity 62 
where acts on ^(Tic) as follows 



L^M ■■= -^Wb, p] + aih - X/e)p{b* - A/6) - - A/e)(6 - A/e), p}. (4.3.3) 



Proof. Using unitary transformation generated by fl4.2.4p of the Hamiltonian, we 
define instead of f l4.1.23|) : 

L^Ap ® Pn) := e*''"i:^,„(e-'^"(p ® pn)e'''")e-'''" (4.3.4) 

= - i[H^, p^pn]+ ae*^"6e-*^"(p ® p„)e^^"6*e-*^" - ^{e'^"b*be-'^" , p ® p„} 
A^ 

= - i[eb*b + {E- — )r/„, p ® Pn] + [(r{b - \/e)p{b* - A/e) 

- |{(^* - ^/e)(& - A/e), p}] ® r/„p„ + {ab*pb - ^{b*b, p}) ® (1 - r/„)p„ 
=L^{p) ® ?7„p„ + L^(p) (g) (1 - ?7„,)p„, , 

here := i^£=o' (14.3.31) . Then discrete evolution operator (14.1.251) gets the form 

C,{p) =Tr^^Je^^-(p® p„)) = Tr«^Je-^^"e^^-(e^^"(p® p„)e-^^")e^^") 

=Tr«^„ (e-^^"e^^-"(e^(''*-^)/Ve~'^''*-'^/^ ® ^nP„, + P ® (1 - r/.)p„)e^^") 

=Tr^^Je-^^"(e-^?(e^(^*-'')Ape-^(^*-'')/^) ® T^^pOe^""") 
+ Tr^,Je-^^"(e-^"(p) ® (1 - r/„)p.)e^^"). 
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Simplifying the last expression 
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+ Tr^^Je"^f(p)®(l-r/„)p„) 

_^g-A(6*-b)AgrLC^^gA(b*-6)A^g-A(6*-b)A^)gA(b*-6)A + (1 _ p)e-^o^ (p) . 

□ 

Now we look for the corresponding adjoint operator C*. It can be calculated using 
the shift transformation (14.3.41) and the adjoint operator f l4.2.12p . For any bounded 
operator B G BiTic) one has: 

£:(i?)=Tr«^Je-^^(i?®p„,)) (4.3.5) 

=pe-A(b-~6)Agr(LC')*^A(6*-6)A^g-A(6*-6)AgA(fe*-b)A _^ (^l - p) e<^0 )* (^B) . 

Here the adjoint operators {L^)* and {Lq)* act as follows: 

(L^)*(B) = z[e6*6, 5] + o-(6* - \/e)B{b - A/e) - |{(6* - A/e)(6 - A/e), B} (4.3.6) 
(L^)*(fi) = (Z^)*(5) = i[eh*h, B] + ^76*56 - ^{6*6, 5} (4.3.7) 

Remark 3. As we indicated in RemarkUl one can extend the adjoint operator C* to 
the algebra of polynomial observables V{b,b*) G 2l('Hc)- 

The proof of the Theorem [8] is similar to the case a = 0. The number of photons 
for the time t = nr can be calculated using the adjoint operator 

N'r^ := a;* (6*6) = Tr^^.(6*6 Cl{pc)) = Tr^,((£:)"(6*6) pc) , (4.3.8) 
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where pc is a gauge-invariant state of the cavity, see Remark |3l 

We look more closely of the adjoint operator £*. Note that if one takes B = b*b 
in fl4.3.5p we get 



Lemma 5. The action of the adjoint operator C* on the number operator of photons 
can be calculated explicitly 

Cl{b*b) =e-""6*6 + p— e^""(l - e''^)b* - p— e-""(l - e''^)b (4.3.10) 
+ p^e"'^^(l-e'^0(l-en- 

Proof Consider the first term in fHX9l) . Let 7A,r(5) = e^(^A)*(5), for B e 2t(Hc). 
To find 7a,t((^* — A/e)(6 — A/e)) we first note that 

(L^)*((6* - X/e){b - A/e)) = iA6 - iXb* - a{b* - \/e){b - A/e) 



and 



, cr, crA , crA 

(L^)*(6) = -266 - -6 + _ = -^6 + _ 



and 



where /i = | + ie. Therefore we have the following system of differential equations 
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for 7A,r 

rf7v((t--A/.)(t-A/.) ^ _ ^^^ji^ _ ^^^^^ ^ ^^^^^^^^j _ ^^^^^^^^.^ 

= -/.7A,.(6) + ^ 

-^r- = -^^'-^^ ^ + 27- 

The solution to this system is 

,,^(,) = - i^) + = e--. + - e-n (4.3.11) 

= - ^) + ^ = + ^(1 - (4.3.12) 

7A ^((f - A/6)(6 - A/e)) = e-"f6 - -6'e-"(-e"' - - + 1) 

e /i /i 

+ A/*-"(^e'=- _ !£ _ 1) + _ e~") + ^e- - ^^^e'^. (4,3,13) 

Making the shift transformation 'jx^riib* — A/e) (6 — A/e)) we get 



e 



-A(b*-fc)A^^^^((5* _ _ A/e))e^(''-'')/^ = e-'^"(6* + A/e)(6 + A/e) 



7 7^: 

A/e(6* + A/e)e-'^"(-e'^^ - - + 1) 
/i /i 



e— -6*6 + — e-""(l - e*'^")6* - ^e-""(l - e-^^fo + —(1 + e""" - 2e-5- coser) 
/i /i 
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Calculating the second term in f l4.3.9p by setting A = 0, we obtain 



+ p^e— (l-e^O(l-en• 



□ 



Note that f l4.3.11|] and (14.3.121) yield 



C*(b*) = e""e-?"6* + p (1 - e'^"e-f n = e-'^'fe* + p— (1 - e~~^^) (4.3.14) 

a — 2ie /X 



and 



(6) = e"*^"e-5"6 - p (1 - e"'^"e"tn = e~^^h - p—il - e"^^). (4.3.15) 

Proof, (of Theorem [8]) The following expression can be found by representing the 
operator £* in the matrix form 

{C*S{h*h) = e-"""6*6 + p— (e"'^"" - e""''")^* - p— (e"""" - e-"'^")6 (4.3.16) 

/i 

+ P^e-'^^(1 - e^O(l - e'n^ p'tTI-T^ (1 - e"^" coser) 

+ P 1 — T^\X ~ ^ ^ COS ner). 



We prove this formula by induction. 
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Suppose (14.3. 16p is true for n, we show it is true for n + 1. 



+ p-^e-^^{l - - en—^ P'fr^V^ (1 " e-^^coser) 



2 

2A2 



+ p —-7,{l-e "2^cosner). 

m 

From the formulas fl4.3.1UI) . 04.3. 14p and 04.3. 15p we get 



{c^Y+\h*h) 

/i /i 

1/^1 

- (e^""" - e""'^")e"''"6 - (1 - e"'^")) 

+ (1 - e''^)(l - e^^)- — -p^-^ . ' _ (1 - e"t-coser) 

1 — e l/ip 1 — e 

+ - e""'^'^cos?T.er). 



1/^ 
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{c*x+\b*b) 

/i ^ 

\2 1 _ -{n+l)aT 

+ p^e-""(l-e'^")(l-e^")- 



-err 



,2A 



2 



'2^^^ ^ (-e-t" cos er + e""" + e-("+^)t- cos(n + l)er - e-("+^)"" 

1 - e-'^^ 



^ I, ,12 1 _ ^-ar 



+ e-("+i)""e-5" cos er - e"("+^)5"e""" cos(n + l)er) 

\2 1 _ -(n+l)ar 

+ p— e-""(l-e'^")(l-e^") 



1-e- 



-P^r-r^ ^^coser) + p^— - e cos(n + l)er) 



which proves (14.3.161) . 

Note that if we take = 1 in (14.3. 16p we get (I4.3.9p . if we take a = we get 
04.2.141) and if r = we get b*b. 

In the limit n goes to infinity we get 

w* - lim {ClY{b*b) 

n—¥oo 

\2 1 9\2 1 9\2 

= p^e— (1-6^(1 -e^^-^ -p^f4.-^(l -e-f^cos6r) 



+ ^2 1 Je-^r i'^ + ^""^(1 - 2P) - 2e— /^(l -p)coser}. (4.3.17) 



□ 

Remark 4. Notice that for a > and p = the limit value ^^^J^ is zero. This is 
trivial because any initial cavity state with finite mean-value of photons ^.S. i| ) will 
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be exhausted by the leaking, a > 0, and the absence of pumping: p = 0. 

Let < p < 1. Then for non-resonant case er ^ 2tts, where s G 1} , the limit 

Yim uc{h*h) = +00 , (4.3.18) 

which corresponds to conclusion of the Theorem about unlimited pumping of the 
perfect cavity, i.e. for a = 0. 

The interpretation of ([^ . 5*. ^ is less transparent in two special cases: 



(a) For the resonant case er = 2txs and a = the mean-value of photons \4-2.1 } is 



bounded and equal to N{0) independent of p. Whereas \4.3.2\j yields the bound 

\2 

YmvuciVh) =p^ ^ , (4.3.19) 

o-->o e 

which is p-dependent. 

(b) For p = 1 the mean-value of photons 1^4. 2.1^ for the ideal cavity is bounded and 
oscillates. Although ( [^ . 5*. ^ for p = 1 gives 

cocib*b) = ^ , (4.3.20) 

for any leaking o" > 0. So clearly, the limits: t ^ oo and a ^ do not commute. 

4.4 Long-time behavior 

Recall that we consider the limiting state (I4.1.28P on the Weyl operator algebra. Our 
result concerning the limiting state of the leaking cavity is the following: 

Theorem 9. Let a > 0. For any gauge-invariant initial cavity state pc and for a 
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homogenous atomic beam with parameter p = Try^j^^{rin Pn), the limiting cavity state 



ucA-) ■= lim 4J-) (4.4.1) 



t~>-oo 



exists and it does not dependent on pc- The explicit form of the limiting functional 



( 4- 1-28 ) on the Weyl operator is 



, ,2 °° 

icr 



ucAW{a)) = e-^ Yl I pexp{- — (1 - e-^^e"^^^"" - ^(1 - e-''^)e'"'^a} + l-p 



k=0 



/i /i 

(4.4.2) 



where := ie + a/2. 



Proof. From Baker-Campbell-Hausdorff formula for the Weyl operators fl2.1.4p and 
from the commutation relations fl2X5|) the action of (L^)* can be expressed as fol- 
lowing 



(L^)*(iy(«)) 

= ie[b*b, W{a)] + a{b* - X/e)W{a){b - A/e) - - X/e){b - A/e), W{a)} 

= {—jjiab — + Jiab* + —(a — a))W^(Q:), (4.4.3) 

^6 



where we denoted /i = | -|- ze. 

Therefore the dynamics generated by (14.3.61) : 

7^^^ := e"(^A)* (4.4.4) 
is quasi-free, which by definition [22] means that for some Tx^c^) and gx^a) 



7A,.(W^(a))=e3^-(")iy(T,,.(a)). 



(4.4.5) 
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See Section [2.2.51 for more details on quasi-free semigroup. 

We can find Tx,r{ci) and gx^rict) using the differential equation the dynamics should 
satisfy 

'^I^^f^ = iLOn,,JW(a))). (4.4.6) 

The right-hand side can be calculated using (14.4.31) equation where instead of a 
we have Tx,t-(«) 



+ ^(T,,,(a) -7VR))W^(TA,.(a)) 

= (-^TA,.(a)6 + 7l7V(")&* -/i|TA,.(«)|' + ^(TA,r(«) 

(4.4.7) 

The derivative of the r-dependent Weyl operator W {Tx^r{.C()) will look like 



2 
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Therefore 7a,t(W^(«)) satisfies the following differential equation 



d7A,r(W^(«)) 



e9A 



c/r c/r ' dr dr 



+ —JZ o^A,r(«) l^+o^A,r(a) 



(ir cir cir 2 ' ' (ir 2 ' cir 

(4.4.8) 

X ^xAW{a)). 

Because of the differential equation for the dynamics (14.4.61) the last equation 
(I4.4.8P should coincides with the equation (14.4. 7p . Then we get the following system 
of differential equations 

dTx,r{a) . . 

1 = -iiTxAoL) 

dr 

and 



dgx^M) 1 .dTx,r{a) 1 . dTx^a) . 2 , ^^/t. ^ ^ ^ / 

= 2^A,r(«) — ^ 2 -^'"^"^ — d^ /^I^A,r(a)| + — {Tx,r{a) -Tx,r{a)) . 

Using the first equation the second one could be simplified and will look like 



-{Tx,ria) - Tx,r{a)) - -|TA,r(«) 



dr 2e' ' " 2 

The solution to the first differential equation is 
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Therefore the second differential equation can be written as follows 

^ = — (e ^^a - e ^^a) e '"^ a r 

dr 2e^ ^ 2 ' ' 

and the solution to this equation is 

= -^(1 - e-") + 1^(1 - e-na - ^(1 - 

Putting the solutions Tx^r{ci) and gx^rij) into the expression for 7A,r(W^(«)) f!4.4.5p 
we get the following 

In order to calculate C*{W{a)) we use f l4.3.5p equation. The first term has the 
following form 



Therefore from f l4.3.5p we get 
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Therefore 



n-l 
fc=0 

(4.4.9) 

In the hmit n goes to infinity W{e~'^'^'^a) converges weakly to 1, so the dependence 
of the hmiting state on the initial state disappears. The characteristic function gets 
the form 



2 _ 

a;c,cr,(Vl/ (a)) = e 2 ||(pe ^ +1— p)VK(e ^ a). 

fe=0 

(4.4.10) 

To see the convergence of the product let us denote 

hk{a) = p[e ^ ^ ^ ' ~ 1)- 

The product Y['k'=o(^~^^i'(^)) converges if and only if the sum J2T=o \^k{c() \ converges. 
Writing hk{a) as a sum of real and imaginary parts we get the following bound 

\ht{a)\ < 2p^(| + 6)(1 + e-f-)e-'=f^ (4.4.11) 

from which the convergence of the series immediately follows. □ 



Remark 5. From ^J^.J^.IO^ and I!i4.4-ll\ ) we see that the state wc,o-(-) is regular (see 



Section l2.1.3\} . since the function a uJc,cr{W{aa)) is C°° -smooth on the vicinity of 
a = 0. By the Araki-Segal theorem (see TheoremU\) , any regular state on CCR('Hc) is 
uniquely defined by the characteristic functional {a i— j- wc,o-(W^(a))}a6C- To check the 
conditions of this theorem we notice that the evolution of the Weyl operator for the 
time t = nr can be written as a convex combination of quasi-free completely positive 
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maps: 

n-l 

m=0 l<ki<...<km<n—l 

X exp (^--^ (/^(l - e-'^")(e-^'^'^" + ... + e-^'"'^")a - c.c.) j . 

Here c.c. stands for a complex conjugate of the first term. From here we see that the 
limiting state ujc.p is an infinite combination of quasi-free states on the Weyl algebra 
and, in general, is not quasi-free itself. 



4.5 Energy flux and entropy production 
4.5.1 Energy variation in perfect cavity 

Since time-dependent interaction in (14.1.21) is piece- wise constant, our system is au- 
tonomous on each interval [{n — l)r, nr). Therefore, there is no variation of energy 
on this interval although it may jump, when a new atom enters into the cavity. Note 
that although the total energy corresponding to the inflnite system f l4.1.3p has no 
sense, its variation is well-deflned. 

Let the n-th atom is actually traveling in the cavity, i.e. t = n{t)T + v{t), with 
n{t) = n — 1 and v{t) G [0, r), see (14.1. lip . Then one can compare the the expectation 
of total energy of the system for the moment t„ = (n — l)r + z/(t„,), with that, when 
the n — 1-th atom was in the cavity, tn-i = {n — 2)r + z/(t„_i). Then by (I4.1.3p . 
(14.1.4P and ( ]4.1.12p the energy variation between two moments tn-i and t„ is deflned 
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by 

AS{tn, tn-l) := Tl{ps{tn)H{tn)) - Tl{ps{tn-l)H{tn-l)) . (4.5.1) 

Lemma 6. For any v G [0, r] one has 

Tr{ps{{n - l)r + u)K^) = Tr{ps{{n - l)r)if„). (4.5.2) 

Proof. 

Tr(p5((n - l)r + u)Hr,) = Tr(e-*'^^"p-s((n - l)r)e*^^"i/„) 
= Tr(p5((n-l)r)i/„). 

□ 

Therefore 

A£'(t„,t„_0 = Tiipsiin - l)r)[iJ„ - H,,_,]) , (4.5.3) 
Let us introduce operator 

7r„,(r) := e*"^"...e'"^V (4.5.4) 

Since fl4.L4p implies 

Hn - H^-i = A (6* + 6) ® {r]n - Vn-i) + l(8)E{r]n- Vn-i) , (4.5.5) 
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by (g^lajl . l^iEM and by [Hk',Vk] = we obtain 

A£{tn,tr,-l) = (4.5.6) 

Tr{pc ® vr„„i(r)(A (6* + 6) ® ]l)<_i(r)[]l ® {r]^ - Vn^i)]} 
+ Tt{pc ® p^(]l ® - Vn-i))} ■ 

Lemma 7. For any n > 1 one gets: 

A " 

7rn{T){b* ® ll)<(r) = e"'^'{b* ® 1) - -(1 - e'"^) J] e^'^"^)'"^! ® r/^ , (4.5.7) 

^ fc=i 
A " 

7Tn{T){b ® Il)<(r) = e-™"^(6 ® 1) - -(1 - e-'^') ^ e'^'^-^)'"! ® 7]k . 

^ fc=i 

Proof. We prove the Lemma by induction. Let 

Bl{r) := e^'^H^* ® ll)e-'"^'= , A; > 1 . (4.5.8) 
Then by fl4.L4p the operator fl4.5.8p is solution of equation 

dsBKs) = i[Hk, Bl{s)] = leBlis) + All ® 77^ , 5*(0) = h*®l, 
which has the following explicit form: 



5*(r) = e'^'{h* ® 11) - -(1 - e'^')l ® rjk . (4.5.9) 

e 



Similarly one obtains 



Bkir) = e-'^\h ® 11) - ^(1 - e-'^')l ® 7]k • (4.5.10) 
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Suppose that fl4.5.7p is true for 7r„, we prove it for vr^+i. 
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A 

k=l 



A 

e ^ 

k=l 

By (14.5.91) and by [Hk',rik] = we obtain 



7r„+i(r)(6* (g)l)= e*("+i)"^(6* ® 1) - -(1 - e^"^)e""^]l ® % 
A 

- -(1 - e'^') V e^'^-i)^"^! ® 
e 

fc=i 



A 



g(n+l)ire^^* (g) 1) - -(1 - e^"') ^ e^*^"^)^""]! ® 



Vk- 



k=l 



Since the similar formula can be obtained for 7r„+i(6 ® the last formula proves 
the lemma. 

□ 

Recall that we suppose that atomic beam is homogeneous, i.e., p = Tr{pc®P^(]l® 
rjn)} is the probability that atom is in its excited state and p is independent of n. 
Then 

Trw^{p^(r/„ir/„J} = Sn,,n2P+ (1 - '^ni.na)^^ ■ (4.5.11) 

Then the second term in the right-hand side of (14.5.61) vanishes. Since we also sup- 
posed that the initial cavity state pc is gauge-invariant, by Lemma [7] and (14. 5. lip 
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one obtains a bound for the flrst term in the right-hand side of f l4.5.6p : 
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Tr{pc ® PA 7r„_i(r)(A {h* + h)® ]l)<_i(r)[]l ® {r]^ - Vn-i)]} 



1 2 



1 - e*^^) J2 e'^'-'^'^'TrnMVn - Vn-i) (4.5.12) 



A 



e 

k=l 

2 '^-1 



-(1 - Yl e-^'~''>'^'TTn^Vk{Vn - Vn-i] 



e 

k=l 



)2 

2 — p(l — p) [cos((n — 2)re) — cos((n — l)re)] 



Hence formulae f l4.5.6p and (14.5. 12p prove for the total-energy variation the following 
statement. 



Theorem 10. The energy variation ( [^ . 5. i ) between two moments t^-i andtn , where 
n > 2 is 



\2 

A£{tn, t„_i) = 2 — p{l-p) [cos((n - 2)re) - cos((n - l)re)] . (4.5.13) 



For the total variation between ti and tn > ti we obtain: 



\2 

A^(t„, ti) = ^^(tk, tk-i) = 2 - p{l-p) [1- cos((n - l)re)] . (4.5.14) 

k=2 



4.5.2 Energy variation in leaking cavity 

In the open cavity the photon energy is not piece- wise constant. Let t„ = (n — l)r + 
uitn) and tn-i = {n — 2)r + z/(t„_i). The energy difference between the times tn-i 
and tn can be divided into three parts: the energy difference between the time 
and the latest time when {n — l)-st atom is still present in the cavity, the energy 
jump when (n — l)-st atom leaves the cavity and n-th atom enters the cavity and the 
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energy difference between the times when the n-th atom is flrst present in the cavity 
and the time tn- 

AE(t„, tn-i) = AE{{n - l)r, + AE{{n - l)r) + AE(t„, (n - l)r). (4.5.15) 

By fl4.1.26p the flrst term can be written as follows 

AE{{n - l)r,t„_i) = TiipsUn - l)r)i7„_i) - Tr(ps(tn-i)i^n-i) 
= Tr(e^^— (Ps((ri - 2)r))i7„_i) - Tr(e'^(*"-)^— (Ps((^ - 2)r))i/„_i). 

Lemma 8. For any < 6 < t 

Tr{e^^'^'-{ps{{n - l)T))Hn) = ee-^^N^dn - l)r) +p— e""^ (4.5.16) 
+ p— e ^-^sme^ + p— ^ 1 -e + p E - — . 

Proof. Using the operator deflned in (14.3.41) the last expression can written as follows 

Tr(e^^-"(p5((n - l)T))Hn) (4.5.17) 

=Tr(p5((^-l)r)e^^-(^n,)) 

=Tr(p5((n - l)r)e-^^"e''^-.-(e*^"i/„e-^^")e^^"). 

From iKT5\i we have that e'^"H„,e-'^" = eb*b + {E - ^)r]n. 

The action of the operator (I4.3.4p can be written explicitly using (14.2.41) . Then 
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continuing the calculation of f l4.5.17p we obtain 

Tr{ps{{n - l)r)e-^(''*-')/^e^(^?)*(e6*6)e^(''*-')/^ ® r/„) 
+ Tr{ps{{n - l)r)e^(^f )*(e6*6) ® (1 - r/„)) 
+ Tr{psiin - l)r)l ® {E - — )r/„) 

= pTtcipciin - l)r)e-^(^*-*)/^e^(^?)* {eh*h)e^^^'-^^'') (4.5.18 



A 



2 



+ (1 - p)ee-'^'^Trc(pc((n - l)r)6*6) + p{E 

by fl4.:ilHp for A = we obtain e^^^o)* (b*b) = e'^^Vh. To calculate the flrst term in 
(14.5. 18p with A 7^ we note that 

{L'iyih^h) = -ah*h+^{h* + h). 



Then we get a differential equation for 7A,r = e 



c?7A,5(fe*fe) _ ^^^^^(5*5) + A^(^,^,(6) + 7,,,(6*)). (4.5.19) 



From (14.3.111) and fl4.3.l2p we flnd 



lx,5{h) + 7v(&*) = e-^'h + e-^V + ^(1 - e"'^') + ^(1 - e"^') 



The the solution of fl4.5.19p yields 
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After the shift transformation, b — j- b + ^, in the last formula we obtain for the 
gauge-invariant initial state pc the A 7^ term 

Trc(pc((ri - l)r)e-^('*-')/^e^(^?)* (er6)e^('*-')/^) (4.5.20) 

\ 2 \ 2 \ 2 

= ee-''^Tic{pc{{n - l)r)b*b) + e-""^— - ^{e''^^ - e"^^) - ^(e"'^'' - e"'^^) 

\ 2 2 

+ ^-^(l + e-'^^-2e-?^cose5). 

Since TTc{pc{{n — l)r)b*b) = N„{{n — l)r) is the number of photons in the cavity 
at the moment t = (n — l)r, by plugging the above expression f l4.5.20p into (14.5. ISp 
we flnd 04.5. 16p . □ 

With the help of this lemma we flnd the energy variance 

AE((n - l)r, = Ti{ps{{n - l)r)i^„„i) - Tr(p5(tn-i)i^n-i) 
= Tr(e^^— (P5((r^ - 2)r))i7„_i) - Tr(e^(*"-)^— (p5((^ - 2)r))i7„_i) 

\2 

= e(e~'"^ - e-'^^(*"-i))iV^((n - 2)r) + p— (e"'"^ - e~'^^(*"-i)) (4.5.21) 

e 

-P7r7^v(e"''" - e-'^^^*"-)) +p^(e-f-siner - e'f sin ez/(t„_i)), 
2e|pr 2 W 

which is the flrst term in the energy difference fl4.5.15p . 
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The third term in (14.5. 15p can be calculated similarly 
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AE{tn, (n - l)r) = Ti{ps{tn)Hn) - Tr(ps((n - l)r)Hn) 
= Tr(e^(*")^-(p5((ri - l)r))i/„) - TiipsUn - 1)t)H^) 
= e(e-'^(*") - i)NJ(n - 1)t) - p—(l - e-''''^'"^) + p^e-^'^'"hmeu(tn) 

(4.5.22) 

To calculate the second term in (14.5.151) note that 

AE{{n - l)r) = Tiipsiin - l)r)if„„i) - Tr(ps((n - l)r)i/„) 
= Tr (e^^— ... e^^-HPc ® PA))(^n - ^^n-i)) 

= Tr (t/:_i(pc ® PA)(^n - i/n-l)) 

= Tr {pc ® PA{K^^y{\{h* + h)® {r]n - Vn-l) + 1 ® ^(?7n - Vn-l))) 

= Tr (pc ® PAiK^iYiKb* + b)(S) (r7„ - r/„_i)) 

= Tr{pc ® p^(t/^„i)*(A(r + 6) ® ]1)[]1 ® ivn - Vn-l)]}, 

where [/^ = f/i'Lnr,s=o = e^'^"'" ••• e^-^"'' was defined in (I4.1.24p . 
Lemma 9. 

{U^Yib ® 11) = e-"^"6 ® 1 - — (1 - e-'^n V e'^'^-^)^"! ® r/^. 

Proof. We prove this lemma by induction. Suppose the formula is true for {U^)* we 
show it is true for {U^_^i)*. The same line of reasoning as in (I4.3.4p can be used to 
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write the action of the operator e"^^^-" 

^ k=l 

By the shift transformation we flnd the right-hand side 



fc=i 



Applying (14.3. lip to the flrst line above we flnd that 

e jj, 
e-*^"+ie"^-"+i(6 ® (1 - r/„+i))e'^"+^ = e-^^^h ® (1 - r/„+i). 

Consequently, 



{K+iTih ® t)) = e-"^^(e-'^^6 - -(1 - e-''^)) ® r^^+i 

A ■ " 

+ e-"'^"e-'^"6 ® (1 - - —(1 - e~^n V e'^'^-^)^^ ® r/^ 

^ k=l 

it 

^ fc=l 



which proves the lemma. 
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Therefore, the energy jump at the moment t = {n — l)r is 

AE{{n - l)r) = Tr{pc ® PA(f/^i)*(A(6* + b)(^l)[l® (r/„ - 



-(1 - e"'^^) ^e"('=-i)^'^Tr(p^77fc(77„ - r/„_i)) 
^ fc=i 



l2- 



The Bernoulh property (14.5.111) imphes 



AE((n- l)r) 

= — (1 - e-'^")p(l - p)e-("-')'^" - —(1 - e-P")p(l - p)e-("-2)'^" 
/i /i 

= 2A2p(l - p)-^i!--^(e-("-2)f- cos - 2)er - e^^"-^)^- cos (n - l)er) (4.5.23) 
+ 4A'p(l - P) ^2^4g2 (^'^""'^^" - 2)er - e'^^-^^^" sin (n - l)er). 



The energy difference between points t„ and is the sum all the three terms 
fl4.5.2m4.5.23D . 

Theorem 11. The energy variation between two moments = (n — 2)r + z/(t„„i) 
and tn = {n — l)r + v{tn), where n >2 is equal to 

AE{tn,t^_,) = e{e-^^ - e^'^'^(*"-))iV<,((n - 2)r) + e(e-'^'^(*") - l)N,{{n - l)r) 

\2 

- p^(l - e""" + e-"^(*"-i) - e-"'^^*")) (4.5.24) 

+ p---7r(e"5^siner + e"^''^*"^ sinez/(t„) - e"^''^*"-!) sin ei/(t„_i)) (4.5.25) 

+ 2A2p(l - p)_ii_(e-("-2)t- cos (n - 2)er - e-(""^)t- cos (n - l)er) (4.5.26) 

(j2 _|_ 

+ 4A'p(l - p) , , (e-("-')t^ sin (n - 2)er - e-("-^)t" sin (n - l)er). (4.5.27) 

(j2 _|_ 



4.5. Energy flux and entropy production 86 
By ( [^ . 5. 24\ ) we obtain for the total variation between ti = and tn > ti : 



AE{tn,ti) =J2^E{tk,tk-l) 

k=2 

n 

= e(e-'^(*") - l)N^{{n - l)r) + e(e-'^^ - 1) ^-((^ - 2 



k=2 

1/^1 1/^1 

\ 2 \ 2 

+ [n — l)p- — —e 2 smer + p-j — rre 2 ^ sm ei/(t„) 
1/^1 1/^1 

+ 2A'p(l - p) , , (e"(""')t^ cos in - 2)er - e-^""^)^" cos - l)er) 
(T^ + 4e^ 

+ 4A2p(l - p) , , (e"(""')t" sin (n - 2)er - e'^^-^^t- sin {n - l)er) 
(j^ + 4e^ 



4.5.3 Entropy production in perfect cavity 

In this section, we make contact with thermodynamics of systems out of equihbrium, 
and in particular with the second law of thermodynamics, see [3]- [5]. 

Let p and po be two normal states on algebra 21(7^). We deflne the relative entropy 
Ent(p|po) of the state p with respect to po by 

Ent(p|po) :=Tr^(plnp-plnpo) >0 , (4.5.28) 

where non-negativity follows from the Jensen inequality: Tr^(pln^) < InTr-^(p^), 
applied to observable A = Pq/ p- 

In the case of the non-leaking cavity we have: 

A5(t) := Ent(p5(t)|ps(t = 0)) = Tr{p5(t)(lnp5(t) - \npc®PA)} , (4.5.29) 



where dynamics is deflned by (14.1.121) . Suppose that all atoms of the beam are in the 
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Gibbs state with temperature which formally can be written as 

PAm:=®PAM. P^-^P) -■= ^Tff^ ^ (4-5-30) 

n>l 

see fl4.1.3p. Since ps{t = 0) = pc®PA = {pc®'^){^<^Pa) and Tr^2®-H^{p5(t) lnp5(t)} = 
Tr-^^,^^^{p5(0) lnp5(0)}, the relative entropy (14.5.291) gets the form 

AS{t): = Tr«,{[pc-P?^]lnpc} (4.5.31) 



13 ^Tr{[p5(0)-p5(nr)](]l®i7^J} 



fe=i 

for t = riT + ly, see (14.1. lip . Here p^*^^ is deflned by fl4.1.15p . 

Remark 6. For any Hamiltonian that acts non-trivially on Tic ® T-Lau ^'^^ ^.i^-U 
Hamiltonian H^^ acting on ^HAk have [Hn,HAf.] = for n ^ k. Note that by 
Iji4-i-14\ ) for our model [Hn, H^f.] = for any n, k. 



Then by virtue of (14.1.150 and (l4.1.16p we have 

Tiipsinr) il(S)HA,)} (4.5.32) 
= TT{e-'^""...e-'^^>'+'ps{kT)e'^^'^+' . . . e'^^" (1 ® Ha,)} 
= TT{ps{ikT) {I^HaJ}. 



By the same argument one obtains also 



Tr{p5(0) {I^HaJ} (4.5.33) 
= Tr{e-*"^^-^..e-'"^V5(0)e'"^i . . . e'"^'^-^ (1 ® Ha,)} 
= Tripsiik - 1)t) (I^Ha,)}. 
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In the case of our model (see Remark [6]) if^^. = e*'^^*(if^^)e~*'^^'''. Therefore the 
last formula gets the form 

Tr{p5(0) {1^H^J} = Tripsikr) (11 ® i/^J. (4.5.34) 

Equations f l4.5.32p and fl4.5.34p shows that the second term in the entropy pro- 
duction fl4.5.3ip vanishes. 

If we suppose that the initial cavity state is Gibbs state for the temperature 1/13 
KTQ^ . then by fl4.1.20p and (14X3|) one gets 

ASit) = TthAIpc - pP] Inpc} (4.5.35) 

= TTuc{[pc-p^c\-(3eb%)} 
= f3e{N{t) - N{0)) , 

where the photon number N{t) is deflned in f l4.2.3p . 

If we denote by AS'"(t) = e{N{t) — N{Q)) the energy variation of the thermal 
cavity deflned by the photon number variation, then f l4.5.35p expresses the 2nd Law 
of Thermodynamics 

AS{t) = (5A£^{t) , (4.5.36) 

for the pumping of cavity. 



Remark 7. In the general situation when [Hn, Hau] combining 31\ ) with 



\4-5.32\j and ( |^.5. 33\) we get for the entropy production at the moment t = nr + v 



AS{t): = TrnA[Pc-pP]\npc} (4.5.37) 

n 

+ /3Y1 Tr{[ps{kT) - ps{{k - l)r)] (]1 ® H^,)} . 



k=l 
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The last term in ^.5.31\) can he rewritten into standard form if one uses the 

identities: 

Tr{ps{kT) (11 ® H^,)} = rr{e^^'= . . . e^^^Pc ® Pa) (1 ® H^,)} 



rr{e-*^^He^^''- . . . e^^^ [pc ® (g) p„,]) ® p,e*^^'= (1 ® i/^J}[]l ® p„,j 

n=l m>fc 



rr{p5((^ - l)r) il®H_^,)}=Trn,^n^^{pt'^^PkHAj • 
For t = fir + i/ie?/ yield the formula for the non-leaking entropy production 

AS{t) = TrnAlPc - pP] In pc} (4.5.38) 

n 

+ Trncm.Mpt'''^Pk)K"A^^H^Je~^-'''' -^^H^.]} ■ 

k=i 
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